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ABSTRACT 

We have undertaken a new ground-based monitoring campaign on the broad- 
line radio galaxy 3C 390.3 to improve the measurement of the size of the broad 
emission-line region and to estimate the black hole mass. Optical spectra and g- 
band images were observed in late 2005 for three months using the 2.4-m telescope 
at MDM Observatory. Integrated emission-line flux variations were measured for 
the hydrogen Balmer lines Ha, H/3, H7, and for the helium line HellA4686, 
as well as <^band fluxes and the optical AGN continuum at A =5100 A. The 
g-band fluxes and the optical AGN continuum vary simultaneously within the 
uncertainties, r cent = (0.2±1.1) days. We find that the emission-line variations are 
delayed with respect to the variable g-band continuum by r(Ha) = 56.3+g; 4 . days, 
t(H(3) = 44.3+|^ days, t(Hj) = 58.1+£? days, and r(Hell4686) = 22.3^ J days. 
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The blue and red peak in the double peaked line profiles, as well as the blue and 
red outer profile wings, vary simultaneously within ±3 days. This provides strong 
support for gravitationally bound orbital motion of the dominant part of the line 
emitting gas. Combining the time delay of the strong Balmer emission lines of 
Ha and H/3 and the separation of the blue and red peak in the broad double- 
peaked profiles in their rms spectra, we determine M^ r = 1.77±g;§? x 1O 8 M and 
using a Hne of the rms spectra M^ r = 2.60±o:3i x 1O 8 M for the central black 
hole of 3C 390.3, respectively. Using the inclination angle of the line emitting 
region which is measured from superluminal motion detected in the radio range, 
accretion disk models to fit the optical double-peaked emission line profiles, and 
X-ray observations, the mass of the black hole amounts to = 0.86TJoig x 
1O 9 M (peak-separation) and = 1.26loi6 x 1O 9 M ((Jn ne ), respectively. This 
result is consistent with the black hole masses indicated by simple accretion disk 
models to describe the observed double-peaked profiles, derived from the stellar 
dynamics of 3C 390.3, and with the AGN radius - luminosity relation. Thus, 
3C 390.3 as a radio-loud AGN with a low Eddington ratio, L eM /L bo i = 0.02, 
follows the same AGN radius - luminosity relation as radio-quiet AGN. 



1. Introduction 

The mass of super-massive black holes (SMBHs) is of fundamental importance for un- 
derstanding active galactic nuclei (AGN) which are powered by mass accretion onto a SMBH. 
Analysis of AGN variability by applying the technique of reverberation mapping (RM) has 
been established as a powerful tool to determine the size of the broad emission-line region 
(BLR) of AGN (Blandford & McKee 1982; Home et al. 2004; Netzer & Peterson 1997; Pe- 
terson 2003) and, under the assumption of virial gas motion, consequently the SMBH mass 
(e.g., Peterson et al.2004). 

Large observational efforts are required to obtain data of sufficient quality and temporal 
coverage to study AGN variability and to employ RM methods. Currently, about 50 AGN 
have been monitored for periods of at least a few months, usually as part of international 
collaborations. The entire available database of reverberation observations obtained through 
2003 was uniformly re-analyzed with improved methods of time series analysis to minimize 
systematic errors (Peterson et al. 2004). Some of these results have been superseded by more 
recent experiments (Bentz et al. 2006b, 2007a; Denney et al. 2006, 2009a; Grier et al. 2008; 
Woo et al. 2010) and additional reverberation results steadily increase the size and quality 
of the database (e.g., Bentz et al. 2009b, 2010a; Denney et al. 2009b, 2010). 
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It has been shown that AGN and non-active galaxies follow the same relation between 
the black hole mass, M^, and stellar velocity dispersion, cr*, of the central spheroidal com- 
ponent of the host galaxy (Ferrarese et al. 2001; Gebhardt et al. 2000; Giiltekin et al. 2009; 
Merritt & Ferrarese 2001; Tremaine et al. 2002) and that the masses based on RM studies are 
consistent with the masses derived from the - a* relation (Nelson et al. 2004; Onken et 
al. 2004; Woo et al. 2010). The close relations of the black hole mass with physical properties 
of the host galaxy, e.g. stellar velocity dispersion, bulge mass and bulge luminosity, indicate 
a coupled growth of the black hole and the formation and evolution of the host galaxy (e.g., 
Silk & Rees 1998; Haiman & Loeb 1998,2001; Bromm & Loeb 2003; DiMatteo et al. 2004; 
Yoo & Miralda-Escude 2004; Volonteri & Rees 2006; Volonteri & Begelman 2010). 

The broad-line radio galaxy (BLRG) 3C 390.3 is a bright and nearby (m v = 15.0, 
z = 0.056; Osterbrock, Koski, & Phillips 1976) FRII radio galaxy with extended double- 
lobed radio emission (Leahy & Perley 1995). After the identification as optical counter- 
part of the radio source 3C 390.3 (Wyndham 1966), it was classified as a N-type galaxy by 
Sandage (1966). Soon, it was discovered that 3C 390.3 shows very broad Balmer emission- 
lines (Lynds 1968) which were later recognized as prominent double-peaked emission-line 
profiles (Burbidge & Burbidge 1971). These double-peaked profiles are generally regarded 
as a characteristic signature of accretion disk emission (e.g., Eracleous & Halpern 1994,2003; 
Gezari, Halpern, & Eracleous 2007). 

So far almost all AGN variability studies have focused on radio-quiet AGN. However, 
3C 390.3 has a well-known variability history (e.g., Cannon, Penston, & Penston 1968; 
Selmes, Tritton, & Wordsworth 1975; Barr et al. 1980; Penston & Perez 1984; Veilleux 
& Zheng 1991; Zheng 1996; Wamsteker et al. 1997; O'Brien et al. 1998; Dietrich et al. 1998; 
Sergeev et al. 2002; Tao et al. 2008), with photometric measurements going back to 1968 
(Yee & Oke 1981). Using the Harvard plate collection Shen, Usher, & Barrett (1972) traced 
back brightness variations to 1895. Thus, 3C 390.3 was a prime target of a multiwavelength 
monitoring campaign in 1994/95 especially because of previous reports of dramatic changes 
in the Balmer line profile shape and strength and to perform coordinated X-ray/UV/optical 
monitoring of a radio-loud AGN for the first time. Currently, 3C 120 and 3C 390.3 are the 
only radio-loud AGN which have been monitored in detail. The multiwavelength study of 
3C 390.3 in 1994/95 (Leighly et al. 1997; O'Brien et al. 1998; Dietrich et al. 1998) shows sig- 
nificant variations in the X-ray, UV, and optical domain over a period of one year. However, 
the light curves of the UV and optical variations show a nearly monotonic increase of the 
continuum and emission-line flux with only moderately strong substructures. The measured 
delays of several broad emission-lines relative to the observed continuum variations are about 
r ~ 20 ± 6 days for H/3 and Ha (Dietrich et al. 1998) and r ~ 36 to 60 days (±18 days) for 
Civ and Lja (O'Brien et al. 1998). More recently, Sergeev et al. (2002) studied the variabil- 
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ity characteristics of 3C 390.3 for a period of nearly one decade. They found a longer delay 
of the response of the H/3 emission with respect to continuum variations compared to the 
result of the 1994/95 study. They concluded that the difference might be caused by different 
continuum variability characteristics which manifest in different continuum auto-correlation 
functions. Recently, Shapovalova et al. (2010), Jovanovic et al. (2010), Popovic et al. (2011), 
and Sergeev et al. (2011) discussed profile variations over periods from a few years up to 
about 15 years. 

In the following work, we present the results of a monitoring campaign undertaken 
at MDM Observatory in late 2005. Based on imaging and spectroscopic data, we find 
that 3C 390.3 clearly showed broad band variations, delayed variations of the broad Balmer 
emission-lines Ha, H/3, and H7, and of the helium line HellA4686. Using the time delays 
and the peak separation in the rms spectra of these emission lines, we estimate a virial 
black hole mass of M^ r = 1.77±oJ? x 1O 8 M for 3C 390.3 and using a Hne of the rms spectra 
M^ r = 2.60±oji x 1O 8 M , respectively. The detection of superluminal motion for 3C 390.3 
indicates an inclination angle of i = 27° ± 2° which allows to correct the observed velocity for 
the tilt of the line emitting region relative to the observer. Taking this into account, we derive 
a black hole mass of M bh = 0.86l£^ x 1O 9 M (peak separation) and M bh = 1.261^;^ x 1O 9 M 
(crune) for 3C 390.3, respectively, which is, within the uncertainties, consistent with recent 
results based on the CallAA8494, 8542, 8662 stellar absorption triplet and the M bh — a* 
relation (Lewis & Eracleous 2006; Nelson et al.2004). 



2. Observations 

We revisited 3C 390.3 to determine the size of the BLR with improved accuracy. In 2005 
from September until December we used the 2.4-m Hiltner telescope at MDM Observatory 
to obtain spectroscopic and imaging data (Table 1). 

The instrumental setup of the Boiler & Chivens CCD Spectrograph (CCDS) and imaging 
camera RetroCam (Morgan et al. 2005) was kept unchanged for the entire observing campaign 
to achieve a homogeneous data set, i.e., avoiding the impact of different instrumental settings. 
The spectrograph was equipped with a Loral 1200x800 pixel CCD, with a projected pixel 
size that corresponds to C.'4/pixel. To cover the entire optical wavelength range, a 150 
grooves/mm grating was used with a fixed grating tilt for the entire observing season. The 
slit width was set to a uniform value of 3" and a position angle of PA= 90°. To allow for 
correction of cosmic ray event contamination, at each epoch two spectra were recorded for 
3C 390.3 with an integration time of tj nt = 900s each. In contrast to previous campaigns we 
observed HD 217086 (BD+61°2373 , 07Vn) for the entire campaign for flux calibration to 
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minimize the uncertainty which is introduced by flux calibrations that are based on different 
standard stars. In addition, HD 217086 was employed to correct for strong atmospheric 
absorption bands, i.e., the 2 C-band (~ 6275 - 6300 A), 2 B-band (~ 6850 - 6950 A), 
and broad H 2 absorption (~ 7160 — 7320 A). This is especially important for the broad Ha 
emission line profile which is severely affected by the B-band absorption. In total, optical 
spectra were recorded for 28 epochs (Table 1). 

In addition, 3C 390.3 was simultaneously observed for 25 epochs in SDSS g-baiad (Fuku- 
gita et al. 1996) using RetroCam at the 2.4-m telescope (Table 1). RetroCam was equipped 
with a E2V CCD 55-20 (1152x770 pixel) with a pixel size of 22.5/mi, which corresponds to 
O'.'259/pixel. The optical design of RetroCam provides the valuable option to switch between 
spectroscopic and imaging mode in less than 2 minutes. Thus, the g-band measurements 
were taken no more than a few minutes before the spectroscopic observations. The exposure 
time was uniformly set to tj nt = 120 sec and generally one exposure was recorded. 

The spectra of 3C 390.3 and the standard star HD 217086, as well as the <^band imaging 
data were processed using MIDAS@ software. Standard dark and flat-field corrections were 
applied to the spectra and g-b&nd imaging data of 3C 390.3 and the HD 217086 spectra. We 
subtracted the night sky intensity for each spectrum individually. The night sky fit was based 
on two regions, ~ 100" to 130" wide on average, which were separated by ~ 60" relative to 
the spectrum of 3C 390.3 and HD 217086, respectively. A 3 rd order polynomial fit was used 
for each wavelength element to fit the spatial intensity distribution of the night sky emission. 

Argon-xenon comparison spectra were recorded for each spectrum of 3C 390.3 for wave- 
length calibration, which is based on ~ 15 to 20 individual lines. The wavelength calibra- 
tion yields a sampling of (3.07 ± 0.05) A/pxl. The spectral resolution, measuring the full 
width at half maximum (FWHM) of isolated lines in the argon -xenon spectra, amounts to 
R ~ 280 at A ~ 5600 A. The spectra of 3C 390.3 were corrected for galactic reddening using 
E B „ V = 0.072 (Schlegel et al. 1998). 

To transform the observed spectra into the rest-frame of 3C 390.3, we measured a redshift 
of z =0.056 ± 0.001 which is based on the narrow emission lines of Ha, H/3, [Ol]A6300, 
[O m]AA4959, 5007, and [Olll]A4363. In the following, we present the analysis of the rest- 
frame spectra of 3C 390.3 and discuss the results based on this spectra. 



Munich Image Data Analysis System, trade-mark of the European Southern Observatory 
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3. Data Analysis 

3.1. (j^Band Imaging 

We used the flat-field corrected images to measure the <^band flux of 3C 390.3 relative 
to comparison stars in the field. The approach of relative photometry has the advantage of 
being quite insensitive to weather conditions, i.e., photometric conditions are not required. 
The BLRG 3C 390.3 is among a sample of 34 AGN which were carefully studied by Penston, 
Penston, & Sandage (1971). They provide Johnson U, B, V magnitudes for three non- variable 
stars in the field of 3C 390.3. RetroCam which has a field of view of 5^0 x 3'.3, allows us to 
observe stars A and B together with 3C 390.3 (Fig. 1). However, with an apparent brightness 
of my = ll. m 71, star A is usually saturated in the observed g-band frames. Star B has a 
magnitude (m v = 14. m 28) comparable to that of 3C 390.3. This star was already successfully 
employed in a prior monitoring campaign of 3C 390.3 as a flux standard (Dietrich et al. 1998). 
Hence, we used star B to determine the relative broad-band flux changes of 3C 390.3. 

To measure the g-band flux of 3C 390.3 and of star B, we applied a square aperture 
with a size of of 13" x 13" (Fig. 1). By using such a large aperture for the comparison star 
and particularly for 3C 390.3 we ensured that the total flux is recorded independently of the 
seeing, which varied from 1"1 to 2"7 with an average of (1"7 ± 0"4) during this monitoring 
campaign. To correct for the night sky flux, the object aperture was separated by a 2" wide 
region (so-called no-man's-land) from a 2" wide region that enclosed the square aperture 
(Fig. 1). This outer region provided an average sky brightness which was used to correct the 
observed flux of the object for the night sky contribution. 

3.2. Optical Spectroscopy 

3.2.1. Intercalibration of the 3 C 390.3 Spectra 

To study broad emission-line flux variations, it is necessary to calibrate the spectra to 
a uniform flux level. This can be achieved using the flux of forbidden narrow emission-lines 
from the narrow-line region (NLR). Due to the large spatial extension of the NLR (of the 
order of 100 pc in Seyfert galaxies) and the long recombination time scales which are due to 
the low gas density of the NLR (r rcc ~ 100 years for n e « 10 3 cm~ 3 ), light travel-time effects 
and long recombination time scales will damp out short time-scale variability. Thus, it can 
be assumed that the flux of narrow emission lines can be treated as constant, at least on 
time scales of years to decades (e.g., Peterson 1993,2000). 
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However, 3C 390.3 needs careful additional attention. The NLR of 3C 390.3 seems to 
be very compact (Baum et al. 1988; Sergeev et al. 2002) and the small emission-line flux 
ratio of F([0 lll]AA4959, 5007) to F([0 m]A4363) indicates the presence of high density gas 
(n e ~ 10 5 cm~ 3 ). Furthermore, narrow-line variability has been reported for this object. 
Clavel & Wamsteker (1987) analyzed the strength of the narrow Lya and CivA1549 emission 
lines. They found that the strength of these two narrow lines decreased by a factor of about 
2 to 3 from 1978 to 1986. On the basis of optical spectra obtained between 1974 and 1990, 
Zheng et al. (1995) presented evidence that the fluxes of [O ill] AA4959, 5007 follow continuum 
variations, although on a longer time scale than the broad emission-lines. During the period 
of decreasing or increasing [O ill] flux there might also be periods of several months or years 
of nearly constant [O ill] flux. Hence, we examined the data closely to test for variability of 
the [O ill] emission lines before we used them for flux calibration. 

To test the reliability of the narrow [O ill] emission lines as an internal calibrator, the 
<^band imaging data are of great value. As noted earlier, these broad band flux measure- 
ments were taken nearly simultaneously with the spectra, i.e., within minutes of the spectral 
observations (Table 1) and the <^band brightness of 3C 390.3 was measured with high preci- 
sion (see §4.2, ~ 0.1% ). To compare the g-band variations to those which can be derived 
from the spectral data, the spectra were convolved with the transmission curve of the SDSS 
<^band filter. The fluxes of the convolved spectra were measured for the transmission wave- 
length range of the g-band filter curve (A = 4000 A to 6000 A). To adjust the flux level of the 
spectra to the g-band measurements, the spectra were rescaled. These scaling factors are on 
average (1.10 ±0.07). 

We used these adjusted spectra, based on the observed g-band flux variations, to test 
whether the [Olll]A5007 emission line flux was constant during this monitoring campaign. 
Although the broad-band flux measurements include the broad and narrow H/3 emission- 
line flux, as well as the [O ill] AA4959, 5007 line emission, the contribution of emission line 
variability to photometric variations can be neglected (Walsh et al. 2009) if it amounts to 
less than ~ 10% of the broad band flux. For 3C 390.3 we find that the H/3 flux amounts 
to less than 5 % of the g-band flux. To test the assumption that the [Olll]A5007 emission 
line flux was constant for the duration of this monitoring campaign, these adjusted spectra 
were internally flux calibrated by scaling each individual spectrum to a uniform flux of 
the [Om]A5007 emission line. First, we calculated a mean spectrum which was used as a 
reference spectrum. Next, we applied a modified version of the scaling algorithm of van 
Groningen & Wanders (1992) to adjust the flux of the [Olll]A5007 in the individual spectra 
to the reference spectrum. This routine minimizes the residuals of the [Olll]A5007 in a 
difference spectrum which is calculated for each spectrum with the reference spectrum. The 
intercalibration to a uniform flux level was achieved by allowing different flux scales, small 
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wavelength shifts, and different spectral resolutions. The average scaling factor amounts to 
(0.99 ± 0.02). This result indicates that the [Olll]A5007 flux was constant in the adjusted 
spectra with a flux of F([0 m]A5007) = (189 ± 4) x fO" 15 ergs" 1 cm" 2 (Table 2). It is 
interesting to note that during the monitoring campaign of 3C 390.3 in f 995 (Dietrich et 
al. 1998) the [Oin]A5007 emission line flux was measured to be F([0 m]A5007) = (161 ± 
4) x 10~ 15 ergs" 1 cm" 2 (rest frame flux) which is about ~ 15% lower. This indicates that 
the [Om]A5007 emission line flux was lower in 1995 compared to 2005 when 3C 390.3 was 
in a less luminous state. Based on the average scaling factor and the small scatter, it is 
justified to assume that the F([0 in]AA4959, 5007) emission-line flux was constant during 
this monitoring campaign. 

While the H/3 - [O ill] AA4959, 5007 region is scaled with respect to the narrow [O ill] A5007 
line, we investigated whether this internal calibration holds for the Ha region at longer 
wavelengths and the fFy - [Om]A4363 region at shorter wavelengths. We intercalibrated 
the spectra for the Ha region with respect to the fluxes of the relatively strong [Ol]A6300 
emission line and the H7 - [Om]A4363 region using the narrow [Om]A4363 emission line. 
We found that the Ha region required small adjustments to achieve a constant [Ol]A6300 
flux. The average scaling factor for the Ha region amounts to (1.014 ± 0.021). For the H7 - 
[Om]A4363 region, minor rescaling of the spectra was also necessary. These scaling factors 
are on average (1.014 ± 0.004). 

3.2.2. Accounting for Host Galaxy and Optical Fell Emission 

A spectrum of an active galactic nucleus is the superposition of various sources, including 
the host galaxy, thermal and non-thermal continuum emission from the AGN, which can 
be described by a power-law continuum, Balmer continuum emission, Fell line emission, 
and emission from other metal line transitions, e.g., carbon, nitrogen, oxygen, magnesium, 
neon, sulfur. These contributions need to be accounted for in order to obtain reliable line 
measurements and the strength of the AGN continuum that most of our study relies on. We 
adopted a multi-component fit approach (e.g., Wills et al. 1985) that we have successfully 
applied in prior studies of high-z quasars and Narrow Line Seyfert 1 (NLS1) galaxies (for 
more details see Dietrich et al. 2002a, 2003b, 2005, 2009). We assume that the spectrum of 
3C 390.3 can be described as a superposition of three components, (i) a power law continuum 
(Fj, ~ u a ), (ii) a host galaxy spectrum (Kinney et al. 1996), and (iii) a pseudo-continuum due 
to merging Fe II emission blends. Neither a Balmer nor a Paschen continuum emission were 
taken into account because the Balmer edge at A = 3645 A is close to the short wavelength 
end of our spectra of 3C 390.3 and the strength of the Paschen continuum is only weakly 
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constrained (Korista & Goad 2001). 

As demonstrated by Bentz et al. (2006a, 2009a) it is crucial to correct for the host 
galaxy contribution to measure the strength of the AGN continuum flux density. Based on 
ACS/HST images which were obtained for 14 AGN of the AGN- watch sample (Peterson 
et al.2004), Bentz et al. (2009a) estimated the host galaxy contribution to the continuum 
flux density at A =5100 A for those AGN, including 3C 390.3. For the large aperture of 
(3" x 9'.'4) which we used, the galaxy contribution amounts to F sa/ (5100A) = (0.87 ±0.09) x 
10~ 15 erg s _1 cm" - 2 A-\ based on ACS/HST imaging data (Bentz priv.comm.). The aver- 
age host galaxy contribution which we determined, based on the multi-component fits, is 
F°^f(5100A) = (0.77 ±0.10) x 10~ 15 ergs -1 cm -2 A -1 , i.e., within the errors consistent with 
the measurements based on the ACS image. To construct representative host galaxy tem- 
plate spectra we combined several individual spectra of various Hubble types. These galaxy 
spectra were retrieved from the publicly available sample published by Kinney et al. (1996). 
The spectral resolution of these galaxy spectra amounts to R ~ 8 A which is comparable to 
the spectral resolution of our spectra. 

To account for Fe II emission, we used the rest-frame optical Fe II template that is based 
on observations of IZwl by Boroson & Green (1992) which covers the wavelength range 
from 4250 A to 7000 A. An alternative optical Fe II emission line template was presented by 
Veron-Cetty et al. (2004) which in addition takes into account NLR contributions to the 
Fell emission. In general, these two templates are similar but in detail they differ at around 
A ~ 5000 A and A > 6400 A. The optical Fe II emission in the spectrum of 3C 390.3 is weak, 
however, and the width of the Fell emission is expected to be quite broad, as is seen for 
other permitted emission lines. Therefore, the choice of the Fell emission template has no 
impact on the results of this study. 

The components (i), (ii), and (iii) were simultaneously fitted to the spectra of 3C 390.3 
that were intercalibrated with respect to the [O ill] AA4959, 5007 emission-line flux, to deter- 
mine the minimum x 2 of the fit. Various host galaxy templates were employed and the best 
results were obtained using an appropriately scaled spectrum of the E0 galaxy NGC 1407, 
which is consistent with the morphological type of the host galaxy of 3C 390.3. The width 
of the Fell emission template was allowed to vary from FWHM = 3000 kms -1 up to 9000 
kms" 1 in steps of AFWHM = 500 kms" 1 . The best fit results were achieved for Fell emis- 
sion templates whose width was on average FWHM = (6000 ± 1000) kms -1 , ranging from 
FWHM = 5000 kms" 1 to FWHM = 8500 kms" 1 . In Fig. 2 we show the fit for a typical 
spectrum of 3C 390.3 from this study, to illustrate the method to determine the different 
components. 

The uncertainty of the fit was estimated based on the distribution of x 2 around the 
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minimum. The spectral slope a of the continuum fit has an average error of Aa ~ 0.01, 
while the uncertainties in the continuum flux density at A = 5100 and the Fell emission flux 
are of the order of ~ 5 % . The best fit of the components (i) to (iii) were subtracted before 
we continued to analyze the emission-line profiles of interest. 

3.3. Narrow-Line Region Emission Lines 

The spectrum of 3C 390.3 shows strong narrow-line emission (Fig. 2). Hence, it is nec- 
essary to correct the broad emission-lines for NLR contributions before line flux and profile 
parameters are measured. In addition, the study of the narrow emission-lines provide infor- 
mation about the ionizing continuum strength. Under the assumption that the NLR emission 
is constant during the observing period in 2005 (September to December), we measured the 
line strength of the narrow emission-lines in the individual spectra of the entire campaign. 
This is an additional test for the quality of the intercalibration of the spectral region of the 
broad emission-lines, studied for variability. 

We used a strong NLR emission line as a template profile to fit the narrow-line spectrum 
and determine its flux (e.g., Whittle 1985c; Dietrich et al. 2005). Using high-quality, high- 
resolution spectra, Whittle (1985a,b) showed that NLR emission-line profiles are similar and 
in particular that the Balmer emission-line profiles are identical to the [O ill] A5007 line profile 
within the uncertainties, even in the case of strong profile asymmetries. Therefore, we used 
the strong emission-line profile of [O ill] A5007 to obtain a representative NLR line profile 
template. To obtain such a profile template we made use of a spectrum of 3C 390.3 which 
we observed with the 2.2m telescope at Calar Alto Observatory /Spain (Aug. 30, 1994). At 
this time 3C 390.3 was at a lower intensity level than in 2005, i.e., the narrow emission-lines 
are more pronounced (Fig. 3). This spectrum which covers a wavelength range comparable 
to this study, has the advantage of a higher spectral resolution (R~ 700 at A ~ 6500 A). 

To isolate the [O ill] emission lines which are located in the extended red wing of the 
double-peaked H/3 line profile we used an appropriately scaled Ha profile to recover the red 
wing of H/3. This scaled Ha profile was subtracted and the residual was used to fit the 
[O ill] AA4959, 5007 line profiles (Fig. 4). The [Om] line profiles are well described using a 
dominant narrow Gaussian component (FWHM = 11.0 A ± 0.1 A) and a broader Gaussian 
component (FWHM = 26.0 A± 1.5 A) which is slightly blue-shifted by v shift ~ 100 kms^ 1 , to 
account for the base of the line profile (Fig. 5). We find that the narrow component carries 
~ 75% of the [Om]A5007 emission-line flux. The individual Gaussian components have a 
priori no physical meaning by themselves, but are just used to obtain a NLR line profile 
template. 
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4. Results 

4.1. The Narrow-Line Region Spectrum 

The [Om]A5007 emission-line profile fit (§3.3) was used to measure the NLR emission- 
line flux of the narrow lines in the spectra of 3C 390.3. While the NLR emission-line profile 
was kept fixed, we allowed the strength of the profile template to vary and also allowed 
variations of the profile width and the location in wavelength space. For measurable emission- 
line doublets, we assumed that the line ratios are F([0 lll]A5007)/F([O ill] A4959) = 2.94, 
F([Oi]A6300)/F([Om]A6364) = 3.05, and F([N ll]A6583)/F([N ll]A6548) = 3.06 (Osterbrock 
& Ferland 2005) in order to fit the weaker line of each doublet. 

Applying this approach we measured the strength of the emission lines and listed them 
in Table 2. The measurement uncertainty for an individual spectrum is of the order of ~ 5 % 
for strong narrow emission-lines like [O ill] AA4959, 5007. In the case of moderately strong 
lines like the [Ol]A6300, [Oll]A3727, [Ne lll]A3869, and blended lines like [N ll]AA6548, 6583 
the error is of the order of ~10% or less, while for the weakest lines it is ~ 25% , e.g., 
[Fex]A6374, [S ll]AA6716, 6731. In Table 2 we present the mean narrow emission-line fluxes 
as measured for all epochs and the scatter of the distribution. 

We corrected the observed line fluxes for internal reddening before the emission-line 
flux ratios were analyzed. Generally, it is assumed that the observed Balmer decrement, 
i.e. F(Ha)/F(H/3) o6s , can be used to estimate the internal reddening (e.g., Davidson & 
Netzer 1979; Netzer 1982). However, there are strong indications that the intrinsic Balmer 
decrement of the emission-line gas of AGN is about F(Ha)/F(H/3) = 3.1 instead of the 
classical case B value for the BLR due to collisional enhancement in the higher density 
partially ionized zone (e.g., Gaskell & Ferland 1984; Halpern & Steiner 1983; Netzer 1982). 

For the NLR gas the strength of the [O n] A3727 emission indicates the presence of a large 
fraction of low density gas (n c p\0 + ) = 3.1 x 10 3 cm~ 3 for 2D 5/2 and nf u (0 + ) = 1.6 x 10 4 
cm -3 for 2D 3 / 2 ; Osterbrock & Ferland 2005). Thus, additional collisional enhancement, par- 
ticularly of Ha emission, is expected to be weak. Therefore, we used the classical value of 
the Balmer decrement of F (H a) / F (H (3) int = 2.85, for pure recombination caseB (e.g., Os- 
terbrock & Ferland 2005) to correct the narrow emission- lines for reddening (e.g., Crenshaw 
et al. 2001). We determine an E B _ V = 0.205 ±0.004 for the NLR of 3C 390.3. Together with 
a Galactic reddening curve (Seaton 1979), we corrected the observed emission line ratios 
relative to F(H/3), applying 

Filine) _ Fjline) 0A Eg _ v {Rx _ 3 _ 68) , . 



- 12 - 



with R\ as the Galactic extinction coefficient at the wavelength of the corresponding emis- 
sion line (Seaton 1979). We also investigated the impact of a larger Balmer decrement, 
F(Ha)/F(H(3) = 3.1 which has been suggested to be more typical for emission line ratios 
of the BLR, i.e., it takes into account collisional excitation and radiative transfer effects - 
(Osterbrock & Ferland 2005). We find that it results in reddening-corrected line ratios that 
are less than ~ 8 % larger for line ratios shortward of H/3 and less ~ 8 % smaller for line ratios 
longward of H/3. The reddening-corrected narrow emission-line measurements (flux cor ), are 
given in Table 2. 

4.2. Continuum Variability 

We measured the relative broad-band flux variations of 3C 390.3 relative to the cali- 
brated comparison starB in the field of 3C 390.3 (Fig. 1). The Johnson V magnitude of this 
star is given as m v = 14.28 (Penston, Penston, & Sandage 1971). During the monitoring 
campaign, starB and 3C 390.3 showed comparable (/-band brightness. This already indicates 
that 3C 390.3 was in a brighter state during this monitoring campaign in late 2005 than in 
1995 (Dietrich et al. 1998). 

Based on the measured <^band flux ratios, we calculate the g-band magnitude variations 
of 3C 390.3. To do this we utilized the known 1^-band magnitude of the comparison starB. 
Fukugita et al. (1996) and Smith et al. (2002) give the transformation of a g' -band into V- 
band magnitude, which involves in addition the (B — V) color. Applying these relations we 
find a (/-band magnitude of m g i = 14.59 (Fukugita et al. 1996) and m g i = 14.62 (Smith et 
al. 2002) for starB, respectively. However, we used the g-band filter of RetroCam (Morgan 
et al. 2005) at the MDM 2.4 m telescope to measure the broad-band flux of 3C 390.3. It has 
been pointed out that there are small differences between the (/band magnitudes (used at 
the 0.5 m telescope of the SDSS) and the g-band filter (used at the 2.5 m telescope of the 
SDSS; Tucker et al. 2006; Davenport et al. 2007; Smith et al.2007). 

Thus, we had to transform the (/-band magnitude of the comparison starB into a 
g-b&nd magnitude. We employed the transformation given by Tucker et al. (2006) with 
g = g' + 0.060 x ((</ — r') — 0.53). To get an estimate of the color g' — r' we used the color of 
the star L107-S97 in Bilir et at. (2008) which has a very similar (B — V) color to comparison 
starB. This results in a correction of A(g-g') = —0.0014, which can be neglected. Therefore, 
for comparison starB the ^band and the <?-,band magnitudes can be treated as identical. 
In Figure 6 (top panel), we show the broad g-band flux density variations of 3C 390.3. 

The errors of the #-band magnitudes of 3C 390.3 and comparison star B are estimates us- 
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ing the uncertainties of the flux measurements. These uncertainties are combined in quadra- 
ture to obtain the error of the flux ratio and hence of the flux density of 3C 390.3. Due to 
the high integrated count rates of comparison starB and of 3C 390.3, in the range of several 
10 6 counts, the errors of the flux ratios are of the order of ~ 0.1 % . 

Next, we measured the continuum flux density variations at A = 5100 A directly from the 
power law continuum fit for each spectrum (Table 3). In Figure 6, the variations of i*\(5100 
A) are displayed. It can be seen that the AGN continuum and the g-band variations follow 
the same pattern and have comparable amplitudes. 

4.3. Broad Emission-Line Flux Variations 

4.3.1. Ha A 6563 

In the spectral region of the broad and double-peaked Ha emission line some of the cor- 
rected spectra still exhibit a weak residual continuum flux level. To adjust for the residual 
continua, we fit a linear local pseudo-continuum. The Ha emission-line flux was measured 
within the wavelength range of 6200 A to 6800 A. This spectral range contains contributions 
of several narrow emission lines, including Ha nar . To correct for narrow line flux contami- 
nation the average intensities of the narrow emission-lines of [O l]AA6300, 6364, [Fex]A6374, 
[Nll]AA6548,6583, [S ll]AA6716, 6731, and Ha nar were subtracted (Table 2) and the resulting 
broad Ha emission line flux is given in Table 3. 

To estimate the uncertainties of the broad Ha flux measurements, first the uncertain- 
ties of the multi-component fit were taken into account. As noted earlier, we assumed that 
the strength of the power law continuum fit, the host galaxy contribution, and the optical 
Fell emission can be determined within 1 %. The next source of error is given by the uncer- 
tainties of the narrow emission-line intensities (Table 2). The largest source of error is the 
level of the pseudo-continuum. The strength of the linear pseudo-continuum fit is based on 
two windows which are located at 6135 A - 6160 A and 6875 A - 6895 A, respectively. Using 
the scatter of the continuum level of these two windows, the error introduced by the pseudo- 
continuum fit was determined. The individual errors were combined in quadrature to obtain 
the uncertainty of the individual broad Ha emission-line flux measurements (Table 3). In 
Figure 6, we present the light curve of the broad Ha emission line flux. Although the scatter 
is larger than for the g-band and F\(5100 A) variations, the broad Ha flux shows a similar 
variation pattern which is shifted to later epochs. 
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4-3.2. H(3 A 4861 and Hell A4686 

As in the case of the Ha emission line, it was necessary to correct for a residual contin- 
uum level of the H/3 A4861 - He IIA4686 wavelength range. The linear local pseudo-continuum 
fits are based on two 20 A wide continuum regions centered at A = 4520 A and A = 5240 A, 
respectively. The corrected spectra were used to measure the broad H/3 A4861 and He IIA4686 
emission-line fluxes. Because the emission-line profiles of 3C 390.3 are so broad (Fig. 2) it is 
necessary to evaluate and to correct for the mutual blending of the H/3 A4861 and HellA4686 
emission, i.e., to identify the contributions of the H/3 A4861 and HellA4686 emission in the 
overlapping region. 

We used the broad Ha A6563 emission-line profile as a template to estimate the shape 
and strength of the broad HellA4686 emission. To extract the broad HaA6563 emission- 
line profile, the mean spectrum of the campaign was used (§4.4) and the mean Ha pro- 
file was corrected for the narrow emission-line contributions ([O l]AA6300, 6364, [Fex]A6374, 
[Nll]AA6548,6583, Ra nar , and [S ll]AA6716, 6731). The width of the double-peaked broad 
emission-line profile of Ha A6563 was scaled to have the same width in velocity space at the 
location of the H/3 A4861 emission line and the strength was rescaled to fit the central part 
and red wing of the H/3 A4861 profile. In Figure 7 we show the H/3 A4861 - He IIA4686 range 
corrected for narrow emission of HellA4686, H/3 A4861, and [O ill] AA4959, 5007. The scaled 
double-peaked Ha emission-line profile describes the broad H/3 A4861 emission profile well, 
except the blue hump which is less prominent in the H/3 profile compared to the Ha A6563 
emission-line profile. This results in an absorption line-like feature in the residual spectrum 
(Fig. 7, bottom panel) which indicates a steeper Balmer decrement for the blue peak in the 
double-peaked profiles. To measure the HellA4686 emission-line flux, the residual emission 
at the location of the He IIA4686 line was fit with a rescaled broad Ha A6563 double-peaked 
profile and also with a single broad Gaussian profile. The profile in the residual spectrum 
of the mean spectrum is better represented by a Gaussian profile than using an appropri- 
ately scaled double-peaked Ha profile (Fig. 7). To determine an appropriate Gaussian profile 
for the broad HellA4686 emission, the blue wing of the Hell line in the residuum was ex- 
tracted and assuming that the broad Hell profile is symmetric the blue wing was mirrored 
to represent the red wing of the HellA4686 line profile (Fig. 8). The resulting profile can be 
described with a single Gaussian profile at A c = 4686.14 A ± 0.04 A and a profile width of 
FWHM = 242.3A±1.4A (15500 km s" 1 ± 100 kms" 1 ). To measure the HellA4686 emission- 
line light curve (Fig. 6) a Gaussian profile of fixed width (FWHM = 15500 kms -1 ) was 
fit to the blue wing of the HellA4686 profile in the spectra of 3C 390.3 (Fig. 8). The flux 
of the HellA4686 emission was measured from this Gaussian fit for a wavelength range of 
A = 4472 A to A = 4900 A. This Gaussian profile fit was subtracted from the spectrum to 
measure the broad H/3 A4861 emission flux in the range from A = 4680 A to A = 5000 A 
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(Table 3). The resulting light curves for the broad H/3A4861 emission line flux, corrected for 
narrow H/3 emission, and the HellA4686 emission line, are displayed in Fig. 6. 

To estimate the total error of these flux measurements, we assumed that the uncer- 
tainties of the multi-component fit are of the order of 1 % . The errors introduced by the 
correction of the narrow emission lines are taken from Table 2. The uncertainty due to the 
broad He IIA4686 emission is of the same order as the total error that is caused by the narrow 
emission lines. As in the case of Ha A6563, the dominant source of error is the level of the 
pseudo-continuum fit. The level of the pseudo continuum fit was varied within the scatter 
given by the two ranges centered at A c = 4520 A and A c = 5240 A that are 20 A wide. The 
individual contributions to the uncertainty were combined in quadrature to obtain the total 
uncertainties of the individual broad H/3 and HellA4686 emission-line flux measurements 
(Table 3). 



4.3.3. H-f\4340 

In the case of the H7 emission line profile region, we correct for a residual continuum 
level as well. The linear local pseudo-continuum fits are based on two 20 A wide continuum 
regions centered at A c = 4160 A and A c = 4530 A, respectively. The corrected spectra were 
used to measure the H7 A4340 emission-line flux (4220 A to 4500 A). These flux measurements 
were corrected for the emission line contributions of the narrow H7 A4340 emission and of the 
[O m]A4363 emission, using the average narrow line fluxes as given in Table 2. The resulting 
light curve of the variable broad H7 A4340 emission line is shown in Fig. 6. The variability 
pattern is very similar to the broad H/3 variations. 

To estimate the total error of the broad H7 flux measurements (Table 3), we assumed 
that the uncertainties of the multi-component fit are of the order of 1 %. The errors intro- 
duced by the correction of the narrow emission lines are taken from Table 2. The dominant 
source of error is the level of the pseudo-continuum fit. The level of the pseudo-continuum 
fit was varied within the scatter given by the two 20 A wide ranges (A c = 4160 A and 
A c = 4530 A). The contributions of the individual error sources were combined in quadrature 
to estimate the total uncertainty of the broad H7 A4340 emission line flux. 
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4.4. Mean and rms Spectra 

44.1. flaA6563, H(3 \A861, and A4340 

Based on the calibrated spectra of 3C 390.3, we calculated mean and residual (hereafter 
simply rms) spectra of the broad emission profiles of Balmer lines Ha, H/3, and H7 which 
are shown in Figure 9. These mean and rms spectra are based on the entire set of MDM 
spectra, except the spectrum that was recorded at the epoch JD = 244 3696. The signal-to- 
noise (S/N) ratio for this spectrum is significantly lower than for all the other spectra of 
3C 390.3 on account of bad weather conditions (S/N2443696 — 10, compared to an average 
signal-to- noise ratio for the sample of S/N = 35 ± 7). 

The mean spectra of the Balmer emission lines show a characteristic double-peaked 
emission-line profile with a strong blue hump while the red hump is very weak. The rms 
spectra of the Balmer lines show a similar shape. The blue hump is the dominant feature 
for the Ha and H/3 rms spectra while it is less pronounced for H7. Only weak residuals of 
narrow emission-lines can be seen in the Balmer rms spectra of Ha, H/3, and H7. The rms 
spectra of the Balmer emission lines also display extended red wings. 

The line width can be characterized by its full width at half maximum (FWHM) and by 
the line dispersion au ne (the second moment of the line profile; e.g., Peterson et al. 2004). In 
order to characterize the uncertainties in each of these parameters we employed a procedure 
described by Peterson et al. (2004). For a set of ./V spectra, iV spectra are randomly selected, 
in particular without regard whether the spectrum has been previously selected or not. These 
./V randomly selected spectra are used to construct mean and rms spectra from which the line 
width measurements are made. This process yields one Monte Carlo realization, and a large 
number (N ~ 10 000) of these realizations yield a mean and standard deviation for each of 
the measurements of the line width. The results of the line profile measurements and the 
corresponding errors are presented in Table 4. 

We find that the mean spectra of the Balmer emission lines show a similar profile width 
of FWHM = 11900 kms -1 to 13200 kms -1 and the second moment of the line profiles <Ju ne 
ranges from 4000 kms -1 to 5400 kms -1 , respectively. The comparison of the rms spectrum 
profile widths show basically the same result. The Balmer emission line rms spectra have an 
average profile width of FWHM = (11430 ± 790) km s^ 1 and an average second moment of 
the line profiles of au ne = (5160 ± 310) kms -1 , respectively. 
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44-2. tfellA4686 

The mean and rms profiles are a Gaussian due to our approach to measure the He IIA4686 
emission line flux and hence also the rms spectrum (§4.3.2). However, the rms spectrum of 
the He IIA4686 variations is contained in the rms spectrum of the broad H/3 emission line 
(Fig. 9). To recover the HellA4686 rms spectrum we examined two approaches. 

First, the fit of the HellA4686 emission-line was subtracted from each spectrum of the 
H/3-[0 lll]AA4959, 5007 emission line region. This results in an uncontaminated broad H/3 
profile. These spectra are used to calculate a mean and rms spectrum which represent the H/3 
variation alone. In Figure 10 the mean spectrum of the H/3-[0 lll]AA4959, 5007 emission-line 
region with and without correction for the broad HellA4686 emission is displayed, as well 
as the comparison of the rms spectra. The difference of these rms spectra can be associated 
with the HellA4686 rms spectrum (Fig. 10). 

As an alternative approach to recover the He IIA4686 rms spectrum, we assumed that the 
rms spectrum of the broad Ha emission-line can be used as a template for the rms spectrum 
of the H/3 emission line. The Ha rms spectrum was re-binned that the width is unchanged 
in velocity space at the location of the H/3 emission-line. In Figure 11 the rms spectra of 
the broad H/3 and Ha emission lines are shown. Next, the re-binned Ha rms spectrum was 
rescaled by dividing by 3.50 to fit the red wing of the H/3 rms spectrum. The difference 
between the H/3 rms spectrum and the scaled, re-binned Ha rms spectrum can be associated 
with the He IIA4686 rms spectrum as well. This difference is displayed in the bottom panel 
of Fig. 11. It can be seen that this representation of the HellA4686 rms spectrum is nearly 
identical to the rms spectrum which is obtained by comparison of the rms spectra of the 
H/3-[0 ill] AA4959, 5007 spectral region, corrected and uncorrected for HellA4686 emission, 
as described above. 

The comparison of the He IIA4686 rms spectrum with the mean Gaussian profile fits to 
the He IIA4686 line profile indicates that the recovered rms spectrum of the broad He IIA4686 
emission line tends to show less variability in the far blue wing than the mean spectrum 
and thus the FWHM is smaller. This result is consistent with the properties of the mean 
and rms spectra of the Balmer emission lines which also show narrower profiles in the rms 
spectra compared with the corresponding mean spectra (Table 4). We used both HellA4686 
rms spectra to measure their FWHM. While the overall shapes of the rms spectra are very 
similar, the signal-to-noise ratio is different. Hence, we used the average of the FWHM, as 
well as the second moment au ne of the profile. The results are given in Table 4. 
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4.5. Variability Characteristics 

We calculated basic variability parameters to characterize the statistical properties of 
the observed emission- line and continuum variations (Fig. 6). The results are presented in 
Table 5. Column (1) gives the spectral feature. The mean flux and the standard deviation 
are presented in column (2) and (3). The excess variance, F var , in column (4) is computed 

as 

F var = ± yfi* - A 2 (2) 

with F as the mean flux (column 2), a F as the rms of the flux variations and A 2 as the mean 
square value of the measurement errors (Rodriguez-Pascual et al. 1997). Finally, in columns 
(5) and (6) the minimum and maximum fluxes are given for each light curve. The relative 
variations of the AGN continua (Fa(5100A) and (/-band) and of the emission lines of Ha, 
H/3, and H7 are of the order of ~ 5 % . The weak He IIA4686 emission line, however, shows 
a relative variation of about ~ 15 % . In spite of the larger uncertainties of the emission-line 
flux measurements of He IIA4686 these uncertainties are accounted for by the parameter A 
in eq. 2 to calculate the relative variability. There is a trend for stronger relative variations 
of emission lines with a higher ionization potential of the involved atom. The larger relative 
variability is also consistent with the results of former variability studies which find that 
F var of high-ionization lines like HeIlA1640 or HeIlA4686 are always larger by a factor of 
about ~ 2 to 3 times than F var for low ionization lines like H/3 or Ha (e.g. Clavel et al. 1991; 
Collier et al. 1998; Dietrich et al. 1993, 1998; Korista et al. 1995; O'Brien et al. 1998; Peterson 
et al. 1991; Reichert et al. 1994; Rodriguez-Pascual et al. 1997; Santos-Lleo et al. 1997; Stirpe 
et al. 1994; Wanders et al. 1997). 



4.6. Time Series Analysis 

Generally, AGN monitoring campaigns have focused on the H/3 - [O ill] AA4959, 5007 
spectral range because the strong forbidden [O m]AA4959, 5007 emission lines can be used to 
establish a robust calibration of the individual spectra. With the large wavelength coverage 
of this campaign for 3C 390.3 we are able to study the delayed response of not only H/3 but 
Ha, H7, and He IIA4686 as well. 

To measure the time delay between AGN continuum variations and variable broad line 
emission, the continuum light curve is used as the driving light curve with the line emission in 
response to these variations. To calculate the time delay between the observed variations, we 
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use the interpolation correlation function (ICCF) method (Gaskell & Peterson 1987; White 
& Peterson 1994; Peterson et al. 1998, 2004). For this monitoring campaign we derived an 
AGN continuum light curve from the power-law continuum fits at A = 5100 A, as well as the 
<^band light curve. 

First, we selected an appropriate continuum light curve. Since the <^band flux contains 
variable broad H/3 emission in addition to the AGN continuum, we studied the impact of the 
broad Hf3 line emission on the g-band variations. The contribution of the host galaxy can 
be assumed to be constant. The g-band light curve and the F A (5100A) light curve (Fig. 6) 
are cross correlated using the ICCF method with F A (5100A) as the driving light curve. The 
resulting CCF is shown in Figure 12. Close inspection of this CCF shows that the g-band 
light curve is delayed with respect to the F A (5100A) variations by r cent = (0.2 ± 1.1) days 
and r peak = (0.6 ± 1.2) days, respectively (Table 6). Within the uncertainties, the F A (5100A) 
variations and the g-band variations are simultaneous and the detected small delay can be 
neglected. It is reasonable to attribute this hint of a delay to the variable broad H/3 emission 
which is part of the g-band flux. Since the scatter in the g-band fluxes is much smaller than 
in the F A (5100A) continuum light curve, we selected the g-band variations as the driving 
continuum in the subsequent cross-correlation analysis. For comparison, we also provide the 
cross-correlation results using Fa(5100A) as the driving continuum light curve (Table 6). The 
cross-correlation functions (ICCFs) for the broad emission lines are displayed in Figure 13. 
The ICCFs show some asymmetric shape and the region around the peak has some structure 
as well. 

To quantify the uncertainties in the time delay measurements, we employ the model- 
independent Monte Carlo FR/RSS method (Peterson et al. 1998), including modifications 
described by Peterson et al. (2004). For each single realization of the method, random subset 
sampling (RSS) is applied in which a light curve with N data points is randomly sampled N 
times without regard to the previous selection of each point. Flux randomization (FR), in 
which a random Gaussian deviation, based on the associated error, is then applied to each of 
the selected N points. This FR/RSS-altered subset of data points is then cross-correlated as 
though it were real data. The peak of the cross-correlation function, r max , which occurs at a 
time lag r pea k, is determined, as is the centroid, r cent . A cross-correlation peak distribution 
(CCPD) for T peak and a cross-correlation centroid distribution (CCCD) for r cent are built 
with a large number (N~ 10,000) of Monte Carlo realizations. We take the average value of 
the CCPD to be r vea k and the average of the CCCD to be r cent . To calculate the uncertainties 
Ar up and Ati ow we assume that 15.87% of the cross-correlation centroid realizations have 
values r > r cent + Ar up and 15.87% have values r < r cent — Ati ow , with analog estimated 
uncertainties for r pea k- These definitions of the uncertainty correspond to ±1 a for a Gaussian 
distribution. 
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4-6.1. Velocity- Dependent Broad Emission- Line Variations 

While the variability characteristics of the integrated broad-line flux can be used to infer 
the size of the BLR, the analysis of the response of different parts of the broad emission- 
line profile, e.g., the line core, and the profile wings, yield information on the dominant gas 
motion, i.e., the kinematics of the broad emission-line gas (e.g., Blandford & McKee 1982; 
Home et al. 2004). Previous investigations of 3C 390.3 have found that the profile wings vary 
in phase (Dietrich et al. 1998; Popovic et al. 2011; Sergeev et al. 2002; Shapovalova et al. 2010) 
which strongly favors orbital motion, i.e., virial motion dominates the gas kinematics. This 
is also supported by accretion disk models to describe the double-peaked Balmer emission 
line profiles (e.g., Eracleous & Halpern 1994; Flohic & Eracleous 2008). Recently, results for 
several AGN have been presented (Denney et al. 2009b; Bentz at al. 2009b) which indicate 
a more complex picture with signatures for orbital motion, as well as for radial motions. In 
addition, for Arpl51 the analysis of spectroscopically resolved broad emission-lines results 
into velocity-delay maps which indicate the presence of substructures like orbiting hot spots 
(Bentz et al. 2010b). 

As in former studies of 3C 390.3, we extracted light curves for different parts of the broad 
emission-line profile. For better comparison with prior studies (e.g., Dietrich et al. 1998; 
Gezari et al. 2007; Shapovalova et al. 2010) the width and placement of those regions are 
motivated by the location and the widths of the blue and red peak in the double-peaked line 
profile. 

In Table 7 the location of the blue peak and red peak for the Ha, H/3, and H7 emission- 
line profiles are given. Within ~ 100 kms -1 both peaks are at similar velocities. Based 
on the mean line profiles (Fig. 9) the width of the peaks amounts to FWHM ~ 50 A and 
FWHM~70A for H/3 and Ha, respectively. We investigated the profiles of the Ha and 
H/3 emission lines only because the H7 emission line is too weak and HellA4686 is best 
represented by a single Gaussian profile. We extracted the light curves of the blue and red 
peak for a 3000 kms -1 wide region. In Figure 14 the boundaries of the blue and red wings, 
the blue and red peaks, and the line center are shown. For comparison the location of those 
profile sections are given in Table 8 which were used by Gezari et al. (2007). The errors for 
the extracted light curves were determined following the same approach as for the integrated 
emission line fluxes. 

Next, the Ha and H/3 line profile section variations are cross-correlated with the g-band 
light curve, as well as with the AGN continuum F A (5100 A). In Figure 15 the corresponding 
ICCF(t) are shown. While for Ha the light curves of the blue wings, the blue and red peak 
and of the line center result into ICCFs with well defined structures, the light curve of the 
red wing of the broad Ha line profile results into a ICCF which provides merely an indication 
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for a possible delay (Fig. 15). For the H/3 line profile the ICCFs are broad and in the case of 
the red wing for K/3 there is only a weak indication for a marginally defined peak (Fig. 15). 

The location of the centroid, the peak and the uncertainties in the ICCFs for the wings, 
the peaks, and the line center of the Ha and H/3 emission lines are determined using the 
technique described above (§4.3) and the results are given in Table 9. As can be seen already 
in Figure 15, the blue and red peak of the Ha line and of the H/3 line vary basically with 
the same time delay relative to the continuum variations, consistent with the model that 
the hydrogen Balmer line emission originates in a rotating disk (e.g., Eracleous & Halpern 
1994; Flohic & Eracleous 2008). Within the uncertainties there is also no delay of the profile 
wings with respect to the peaks or the center of the line profile. To test this result we also 
applied the time series analysis to the variations of the blue and red wings and the blue and 
red peaks to eliminate uncertainties which are introduced by the continuum light curves. 
The corresponding ICCFs for the Ha and H/3 lines are displayed in Figures 16. Within the 
uncertainties we detect no time delay between the variations of the blue and red peaks and 
wings (Table 9). 



4.7. Stochastic Process Estimation for AGN Reverberation 

Generally, the time delay between continuum and broad emission line flux variations 
is determined by calculating cross-correlation functions. However, for a better specification 
of the possible impact of gaps in the temporal coverage of the light curves, recently an 
alternative approach has been suggested by Zu et al. (2011) and applied by Grier et al. (2012) 
to the measured variability of the NLS1 galaxy Mkn335. This method to analyze irregularly 
sampled light curves with a more statistical approach was developed by Press et al. (1992) and 
Rybicki & Press (1992). They assumed that an irregularly sampled light curve can be treated 
as a damped random walk process with a specific amplitude a and an exponential damping 
time scale r^. In recent years it has been shown that quasar light curves can be, indeed, 
well described with this approach (Kelly et al. 2009; Zu et al. 2012) and furthermore that 
the amplitude a and damping time scale Td for quasar variability are found in a well defined 
region of the a - Td space (Kozlowski et al. 2010; MacLeod et al. 2010). Recently, Mushotzky 
et al. (2011), however, have shown using Kepler data for four AGN, that the variability of 
all AGN cannot be described with the random walk approach but that individual AGN have 
intrinsically different variability characteristics. 

To investigate also AGN variability data and in particular the impact of gaps in the tem- 
poral sampling in AGN light curves Zu et al. (2011) have modified this statistical approach. 
We used their method of Stochastic Process Estimation for AGN Reverberation (SPEAR) 
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as presented by Zu et al. (2011) to investigate the variations of 3C 390.3 (Fig. 6). First, we 
used the g-band variation and the F A (5100A) variability to test whether the SPEAR results 
are consistent with our cross-correlation results (Table 6). We started with an amplitude of 
a = 0.05 and an damping time scale of = 60 days. With the g-band variations as a driver 
light curve (because of its much smaller errors than the F A (5100A) light curve) , we find using 
SPEAR a delay of r = 0. 55^45 days which is consistent with the cross-correlation result. 

Next, we used the g-band light curve as the driver to recover the time delays of the 
Ha, H/3, H7, and HellA4686 variations. Using SPEAR we computed time delays for the 
strong Balmer lines Ha and H/3 which are basically identical within the errors to the cross- 
correlation based delays. In the case of the H7 and HellA4686 emission lines, however, 
the delays determined with SPEAR and based on the cross-correlation analysis show some 
differences. The time delays of the Balmer lines Ha, H/3, and H7 derived with SPEAR 
are consistent with those of the ICCF analysis using F A (5100A) as the driving continuum. 
However, for the He IIA4686 emission lines the delays are just consistent within the l-a errors. 
The time delay of the H7 variations based on the g-band light curve as the driver continuum 
and the SPEAR result differ by nearly a factor of 2 (Table 6). Furthermore, due to the larger 
uncertainties of the flux measurements for the H7 and HellA4686 emission lines, SPEAR 
needed to be restricted to a narrow range of allowed time delays to avoid spurious delays at 
t ~ 120 days, which is longer than the period covered by the experiment. 

SPEAR treats large gaps in the temporal coverage of a light curve in a statistically 
based approach and also follows the impact on the uncertainty of the time delay. However, 
in the case of the study of 3C 390.3 with no large gaps in the light curve, we applied SPEAR 
predominantly to compare the results with those of the ICCF analysis and the results of both 
methods do well agree. Hence, we used the results of the ICCF analysis in the following 
for consistency with former variability investigations which allows an easier comparison with 
those results. 



4.8. Black Hole Mass and Eddington Ratio 

4-8.1. Black Hole Mass Estimates 

To determine the mass of the black hole in the center of 3C 390.3, we assumed that the 
line-emitting gas is bound, i.e., that the virial theorem can be applied. This approach is well- 
motivated by the results of the fitting of the double-peaked emission line profiles of 3C 390.3 
with an elliptical disk model (Flohic & Eracleous 2008). In addition, for several AGN it has 
been possible to measure the time delay of additional broad emission line flux variations. 



These studies consistently show that emission-lines with higher ionization energies show 
shorter time delays, in keeping with the observed ionization stratification of the BLR and 
their closer locations to the central continuum source (Bentz et al. 2009a; Kollatschny 2003; 
Onken & Peterson 2002; Peterson & Wandel 1999,2000). Together with broader emission 
line profile widths, these results provide strong support that the gas is in gravitational bound 
motion. The virial mass of the black hole is given by 



with cr as the size of the broad emission line region, Av which is associated with the line 
profile width, and G as the gravitational constant. The factor / depends of the geometry, 
kinematics, and the orientation of the line emitting region. It has been shown by Onken et 
al. (2004) that the average of / amounts to < / >= 5.5 when the virial black hole mass, 
based on a rms spectrum profile width, is calibrated to the M^h — a* relation for quiescent 
galaxies (e.g., Tremaine et al. 2002; Gultekin et al. 2009, 2011). In the following we calculate 
the pure virial product, i.e., we will assume / = 1.0 because in the case of 3C 390.3 and 
the successful accretion disk models v Q i> s = Vi ntr x sin? can be used to derive the intrinsic 
velocity, v intr , of the gas. 



Using the second moment of the emission-line profiles au ne (cf. Peterson et al. 2004) for 
the mean and rms spectra (Table 4) and the measured time delay of the broad emission-line 
flux response to continuum variations, given by r cen t (Table 6), derived from a threshold of 
80% of the maximum of the ICCF, we calculated the virial black hole mass using eq.(4). 
The results are presented in Table 10. For comparison we estimated black hole masses for 
the individual emission lines using au ne from the the mean and the rms spectra and the time 
delay which we derived using the AGN continuum variations of F A (5100A) and of the g-hand, 
as well (Table 6). We find that within the uncertainties the virial mass, M^ r , of the black 
hole of 3C 390.3, based on the mean emission line profiles, is in the range of 1.0 x 1O 8 M 
to 2.6 x 1O 8 M based on the properties of the mean profiles of the Balmer lines and the 
HellA4686 emission line, independent of the continuum used (Fig. 17, Table 10). However, if 
only the strong Balmer emission lines Ha and H/3 are considered because the uncertainties for 
the H7 and HellA4686 emission lines are larger, the virial products are in better agreement, 
2.3 x 1O 8 M to 2.6 x 10 8 M Q . The rms spectra yield results which are consistent with the 
mean-spectrum with virial black hole masses of 1.1 x 1O 8 M to 3.0 x 1O 8 M (Table 10) 



M bh = f 



ct Av 2 



(3) 



G 
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and the strong Balmer emission lines Ha and H/3 yield black hole mass estimates of about 
~ 2.6 x 1O 8 M . Since, the rms spectrum represents the actually variable part of the emission 
line profile, we conclude that the virial product of the black hole of 3C 390.3, based on the 
strong Balmer emission lines Ha and H/3 is M^ r (3C390.3) = (2.60t°;^) x 1O 8 M . 

To estimate the black hole mass using reverberation mapping studies, it is still nec- 
essary to specify the factor / in eq. (3) which depends on the not well-constrained geometry 
and kinematics of the gas, as well as on the orientation of the line emitting region. While it 
is generally assumed that this factor / is of the order of unity, it has been pointed out that 
a small inclination of the line emitting region can result in a significant underestimation of 
the black hole mass (e.g., Collin et al. 2006) if the gas motions are almost entirely in the disk 
plane with little motion in the polar direction. In contrast to almost all AGN monitored so 
far, the BLRG 3C 390.3 is unique since it is the only radio loud AGN with extended double- 
lobed radio emission (Leahy & Perley 1995) whose variability has been studied in detail. 
This allows us to estimate the impact of the inclination % of the line emitting region on the 
black hole mass measurement and hence the factor /, assuming pure disk motion of the gas. 
Based on VLBI observations, Alef et al. (1988) reported indications for super-luminal motion 
for 3C 390.3. Further radio observations and a detailed analysis provided information on the 
inclination i of the axis of the radio emission in 3C 390.3 (Alef et al. 1994, 1998), which was 
found to be i = 28° +go. The inclination of the radio axis of 3C 390.3 has been also estimated 
based on X-ray observations (Eracleous et al. 1996). They successfully modeled the resolved 
Fe Ka line emission with X-ray reprocessing in a cool, dense disk of gas at an inclination 
of i ~ 26°. A third way to estimated the inclination of the disk-like line-emitting region is 
provided by fitting the double-peaked profile. Using a relativistic Keplerian disk Eracleous 
& Halpern (1994) derived an inclination angle of i — 26° ^2°- Recently, Flohic & Eracleous 
(2008) analyzed the variations of the double-peaked profile shape of Ha for Arp 102B and 
3C 390.3. They found that the best disk models for 3C 390.3 strongly indicate an inclination 
angle of i = 27° ± 2°. 

In a simple model of an accretion disk as the dominant source of the observed double- 
peaked emission-line profile, the observed velocity is actually v i> s = Vi ntr x sini Assum- 
ing an inclination angle of i — 27° ± 2° this implies a correction factor of / = 4.85lo'gQ. 
Within the uncertainties, this factor is consistent with < / >= 5.5, determined by Onken 
et al. (2004) and < / >= 5.25 by Giiltekin et al. (2009) which were determined by compar- 
ing the reverberation-based black hole masses with the masses determined using the M-a 
relation, with the stellar velocity dispersion a based on stellar or gas dynamics of the host 
galaxy's bulge. Applying the correction due to the measured inclination of the accretion disk, 
the black hole mass of 3C 390.3 based on the Ha and H/3 rms spectrum profile properties 
and r cent yields M bh (3C 390.3) = (1.26+°;^) x 1O 9 M . 
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Accretion disk models have been shown to be successful in explaining the observed 
profiles and even profile variations over longer time scales for 3C 390.3, Thus, the typical 
velocity of the line emitting gas can be associated with the location of the blue and red 
peak of the double-peaked emission-line profile, i.e., the location in velocity space of the 
dominant part of the line emitting gas. Using the separation of the blue and red peak 
of the Balmer emission-line profiles as Av, we estimated the virial product using eq.(4) 
and r cent (80% threshold level). The derived black hole masses are given in Table 10. 
Based on the Balmer emission lines Ha, H/3, and H7 we find virial black hole masses in 
the range of 1.1 x 10 8 £ M bh £ 2.0 x 1O 8 M using the mean emission-line profiles and 0.9 x 
10 8 < M bh < 1.9 x 1O 8 M based on the rms spectra (Fig. 17, Table 10). However, this wide 
range is caused by the H7 based results which have a larger uncertainty than the black hole 
mass estimates obtained from the strong Ha and H/3 emission lines which are consistent 
within the errors. Since the rms spectrum represents the actual variable part of the emission 
line gas, we calculated the average virial black hole mass which is based on the rms spectra 
of the Balmer emission lines Ha, H/3, and H7. We find a virial black hole mass of M^ r = 
(1.77lo;3i) x 1O 8 M . Together with the correction factor / = 4.85±qH which is given by 
the inclination of the accretion disk, we find for the black hole mass for 3C 390.3 a value of 
M bh = 0Mt° fs x 10 9 M . 



4-8.2. Eddington Ratio 



Using the black hole mass estimates based on the broad Balmer and emission-line profiles 
in the rms spectrum and the optical continuum luminosity of 3C 390.3 at A = 5100 A, we 
computed the Eddington ratio L bo i/L edd . The value of the conversion factor f L between the 
monochromatic luminosity and the bolometric luminosity (L bo i — f L x AF A (5100)), is still 
debated in the literature (e.g., Elvis et al. 1994; Laor 2000; Netzer 2003; Richards et al. 2006). 
Recently, Marconi et al. (2008) suggested even a luminosity-dependent correction factor fi, 
due to radiation pressure effects. We assumed that the bolometric luminosity is given by 
L bol = 9.74xAL A (5100A) (Vestergaard 2004). We find = 2.22x 10 45 ergs" 1 . To calculate 
the Eddington luminosity L edd , we assumed that the gas is a mixture of hydrogen and helium 
(fj, = 1.15), i.e., L edd = 1.45 x 10 38 M bh /M ergs -1 . The derived Eddington ratio L bo i/L edd 
for 3C 390.3, based on this monitoring campaign, amounts to L bo i/L edd = 0.0181q;qo5- This 
low Eddington ratio is typical for radio-loud AGN like 3C 390.3 (e.g., Boroson & Green 1992; 
Boroson 2002). 
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5. Discussion 

5.1. Times Series Analysis and the Size of the BLR 

Recently, results of studies on the long-term variability properties of 3C 390.3 have 
been presented by Sergeev et al. (2002,2011) and Shapovalova et al. (2010). In two studies 
Sergeev et al. (2002, 2011) investigated the correlated variations of the optical continuum 
and the response of the broad H/3 emission-line flux for the years 1992 to 2000 and 2000 
to 2007, respectively. In both studies they found that the H/3 variations are delayed by 
t(H/3) = 82+^ days (1992 - 2000) and r(Hf3) = 94 ± 6 days (2000 - 2007) and for Ha 
delays of r(Ha) = 162±?| days (1992 - 2000) and r(Ha) = 174 ± 16 days (2000 - 2007) were 
determined. A comparable result was found by Shapovalova et al. (2010) who studied the 
optical variations for 3C 390.3 from 1995 to 2007. They report that the variations of the broad 
H/3 and Ha emission-line flux are delayed by r(H/3) = days and t{Hq) = 127±18 days, 
both too long to be measurable from our data. However, due to the lower sampling rate of 
the Ha light curve the ICCF analysis displays two possible peaks also at r(Ha) ~ 24 days 
and r(Ha) ~ 151 days, respectively. 

Those results are inconsistent with the delays which were reported by Dietrich et 
al. (1998) using the data of the X-ray- UV- optical monitoring campaign in 1994/95 (Leighly 
et al. 1997; O'Brien et al. 1998; Dietrich et al. 1998), i.e., t(H/3) = 23 ±4 days and r(Ha) = 
19 ± 9 days and also with the results of this study (Table 6). However, the time delays 
which we determined for 3C 390.3 in 1995 are consistent with the delays we find for Ha and 
H/9 in this study. In 1995 the strength of the optical continuum corrected for host-galaxy 
contributions was Fa(5100A) = 1.16 x 10" 15 ergs -1 cm -2 A" 1 . During the 2005 monitoring 
campaign the continuum was about ~ 6x stronger (Table 5). Using the radius - luminosity 
relation for AGN (Bentz et al. 2009a) it can be expected that the derived delays of the Ha 
and H/3 emission-lines should be about ~ 2.5 x as long than those of 1995. This is in good 
agreement with the delays which we have measured (Table 6). 

To investigate the cause of the discrepancy between the results based on the analysis of 
the variations over about 10 years compared with those covering about three months up to 
one year, we re-analyzed the data published by Sergeev et al. (2002,2011) and Shapovalova 
et al. (2010). Using the data of the Sergeev et al. study, we compiled H/3 and F A (5100A) 
continuum light curves which cover about 15 years with 413 epochs for the continuum and 
131 epochs for H/3. We applied the same ICCF analysis to these light curves as we did 
for the measurements for this study. We found r(Hf3) = 8I.6I22 5 days which is consistent 
with the results of Sergeev et al. and due to the large uncertainties, within 1-a it is even 
in the range with the result of this study. Next, we studied parts of the light curve which 
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covered approximately 1000 days to study the impact of the luminosity state of 3C 390.3 on 
the delay. We found delays of the broad H/3 variations ranging from t{H(3) = 13+^ days 
to t(H(3) = 102^}gdays with no clear relation between the continuum state and the delay. 
However, the mean spacing, especially for the H/3 emission-line flux measurements, is only 
~ 20 days to ~ 160 days, i.e., probably not sufficient. Only for the last three years of the 
light curve the temporal sampling is about ~6 days for the continuum and ~20 days for the 
H/3 emission-line flux. For this period we find a delay of t(H/3) = 102+gdays. The H/3 
emission line flux measurements by Sergeev et al. (2011) and by Shapovalova et al. (2010) 
do overlap with our monitoring campaign. Therefore, we can compare those with our H/3 
emission line light curve (Fig. 18). It can be seen that the continuum light curve follows the 
variations we measured while the coverage in H/3 is not as high. 

In the same way, we re-analyzed the continuum and H/3 emission-line light curves given 
by Shapovalova et al. (2010). For the entire data set which covers about 12.5 years we 
find r(Hj3) = 96. 2tH'% days, consistent with the results of Sergeev et al. Similar to the 
Sergeev et al. data the light curves were split up into periods of about 1000 days duration 
with temporal sampling of about 40 days. We found delays of the broad H/3 variations 
ranging from r(H(3) = 16^ days to r(Hf3) = 210^40 days with no clear relation between 
the continuum state and the delay. 

In a study on the reliability of cross-correlation function time delay determinations 
Welsh (1999) has pointed out that variations on longer time scales, e.g., on the dynamical 
time scale of the BLR which is of the order of years, will shift the response time to longer 
delays. Long term variations are different from direct reverberation signals which measure 
the instantaneous response of the emission-line gas to continuum variations, while long term 
variations trace a more gradual response to an overall increase or decrease of the contin- 
uum strength. To correct for this effect it is necessary to detrend light curves from those 
gradual changes. To detrend the continuum and H/3 emission-line variations, we fit a low 
order polynomial to the light curves. Next, we applied the analog time series analysis to 
the detrended light curves of the Sergeev et al. (2002,2011) and Shapovalova et al. (2010) 
studies. We derived time delays of the H/3 variations of r(H(3) = 90.4^g2days (Sergeev et 
al.) and r(Hf3) = 82.7^}|g days (Shapovalova et al.). Within the uncertainties (1-a errors) 
these delays are consistent with those obtained with the original light curves. However, the 
uncertainties of the time delays using the detrented light curves are about a factor 2 to 3 
smaller than those for the light curves including gradual long term variations. 

Furthermore, the properties of continuum strength variations have a significant effect 
on the time delay measured from cross-correlation functions. This issue can be responsible 
for the different delays measured in our study and those using variations of the continuum 
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and the broad emission lines over more than 10 years, in particular, the auto-correlation 
function of the continuum, ACF cont . It has been already noted by Sergeev et al. (2002,2011) 
that the width of the ACF cont of the continuum light curve which is covering several years 
up to more than a decade is much broader than the ACF cont of a shorter campaign like the 
one in 1994/95 and that this will result in a longer time delay (see also the Appendix for a 
more detailed discussion). 

In addition, our campaign covered only a little more than 80 days and hence we are 
not able to measure time delays of about 100 days and more. We think that the different 
time delays are caused by the long time period of about a decade to cover the variation of 
3C 390.3 and the wider temporal sampling of the measurements. Furthermore, the signifi- 
cantly different widths of the ACF cont and the fact that in addition to reverberation signals, 
i.e., the direct response of line emitting gas to continuum variations, also variations of the 
emission-line flux are included. These are associated with changes of the physical conditions 
and the distribution of the gas which happen on dynamical time scales and which are uncor- 
rected with continuum variability also contributes to the different time delays. Therefore, a 
long duration campaign for 3C 390.3 with densely sampled measurements will be necessary 
to find a definitive result. 



5.2. Black Hole Mass Estimates 

A wide range of mass estimates for the super-massive black hole of 3C 390.3 have been 
reported, ranging from ~ 1.3 x 10 8 M up to ~ 7 x 10 9 M using emission line profile proper- 
ties and estimates of the size of the BLR (e.g., Barr et al. 1980; Bentz et al. 2009b; Clavel & 
Wamsteker 1987; Gaskell 1996; Peterson et al. 2004; Sergeev et al. 2002,2011; Wamsteker et 
al. 1997), as well as using the Call triplet in the near infrared employing the M-cr* relation 
yielding about M bh ~ 4 to 5x 10 8 M (Nelson et al. 2004). Most of these studies favor a 
black hole mass for 3C 390.3 of the order of 5 to 10 x 10 8 M Q . Our measurement of the mass 
of the black hole of 3C 390.3 with M bh = 0.86+^^ x 1O 9 M (based on the separation of 
the blue and red peak in the rms spectrum) and M bh = 1.261°;^ x 1O 9 M (using ou ne ) 
is consistent with the black hole mass based on the M — a* relation. Using the M-a* re- 
lation from Giiltekin et al. (2009), we calculated the stellar velocity dispersion which is 
expected for the host galaxy of 3C 390.3 based on our estimated black hole mass. We find 
er* = 257^ 9 Qkms~ 1 for 3C 390.3 which is consistent with <r* = 273 ± 16kms _1 as measured 
by Nelson et al. (2004). Furthermore, simple disk models which we applied to describe the 
overall structure of the broad double-peaked hydrogen Balmer emission lines yield a mass 
estimate of M bh ~ 10 9 M . 
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An additional test for the reliability of the derived black hole mass for 3C 390.3 is given 
by the comparison of our measured values for the size and continuum luminosity of 3C 390.3 
with the expected values derived from the radius -luminosity relation. Guided by simple 
photoionization models, a relation between the continuum luminosity of an AGN and the 
radius of the BLR is expected and Kaspi et al. (2000) provided the first convincing evidence 
for such a relation. Careful re-analysis and additional observations (e.g., Peterson et al. 2004; 
Bentz et al. 2006b,2007a,2009b; Denney et al. 2006,2009a,2010; Grier et al. 2008; Onken et 
al. 2003) have reduced the uncertainty of the slope of the relation. It turned out that the 
correction for host galaxy contamination has a profound impact on the R L relation (Bentz 
et al. 2006a,2009a). With a slope of a = 0.52 ± 0.04 (Bentz et al. 2009a), we estimated the 
black hole mass of 3C 390.3. For the continuum flux at A = 5100 A we used the average 
continuum flux of the AGN continuum (Table 5), corrected for host galaxy contributions, 
with AL^(5100A) = 2.28 x 10 44 ergs" 1 and the measured delay for the H/3 emission line is 
T cent — 44 days. Using the radius - luminosity relation as given in Bentz et al. (2009a) 



log Rblr = K + a\og(X F A (5100A)) (5) 

with K = -21.3l|i and « = 0.519laojjjj the H/3 BLR radius amounts to r = 52.7±|;1 
days, which given the intrinsic scatter in the relationship, is consistent with the measured 
r cent (Hf3) = 44.311;° days ™ the r cent (H/3) = 47.91;;; 4 days using SPEAR. 



6. Summary 

We present results of a ground-based monitoring campaign on the broad-line radio 
galaxy 3C 390.3. Optical spectra and <^band imaging were obtained in late 2005 for three 
months using the 2.4-m telescope at MDM Observatory. Integrated emission-line flux varia- 
tions were measured for the Balmer lines Ha, H/3, H7, and for the helium line HellA4686, as 
well as g-hand fluxes and the optical AGN continuum at 5100 A. The g-band fluxes and the 
optical AGN continuum are varying simultaneously within the uncertainties (r cent = —0.2 ± 
1.1 days). We measure time delays for the emission-line variations with respect to the variable 
#-band continuum of r(Ha) = 56.3l£|days, r(H(5) = 44.3l^days, r(H-f) = 58. ll^ days, 
and r(Hell4686) = 22.3l3g days. The Dme an d red peak in the double peaked line profiles, 
as well as the blue and red outer profile wings vary simultaneously within ±3 days. This 
provides strong support for gravitationally bound orbital motion for the dominant part of 
the line emitting gas. Using the separation of the blue and red peak in the broad double- 
peaked profiles in the rms spectra of the Balmer emission lines and the corresponding time 
delays we determine a virial black hole mass of M^ r = 1.77^^1 x 1O 8 M for the black hole 
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of 3C 390.3. Using the inclination angle % = 27° ± 2° of the line emitting region the intrinsic 
velocity, Vi ntr , can be recovered from the measured v obs = v i ntr sin %. This results in a black 
hole mass of M bh = (0.86^) x 1O 9 M for 3C 390.3 and M bh = (1.26±g;?J) x 1O 9 M based 
on <Ju ne of the rms-spectrum. This mass estimate is consistent with the mass indicated by 
simple accretion disk models to describe the observed double-peaked profiles, as well as with 
black hole masses derived from studies on the stellar dynamics of 3C 390.3. Furthermore, the 
mean continuum luminosity and the measured time delay for the broad emission-line flux 
variations of H/3 is consistent with the most recent AGN radius - luminosity relation. Thus 
3C 390.3 as a radio-loud AGN with a low Eddington ratio of only L edd /L bo i = 0.02 follows 
the same AGN radius - luminosity relation as radio-quiet AGN. 

We thank J. Halpern, S. Tyagi, and all the observers at MDM Observatory who con- 
ducted the observations in Fall 2005, for the first time in service mode at MDM. We also 
acknowledge financial support from NSF grants AST-0604066 and AST-1008882 to OSU. 
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A. Auto-Correlation Function and Time Delay 

Under the assumption that the observed emission-line flux L(t) is the superposition of 
the response of the emission-line gas to continuum variations C(t) of gas with the same 
time-delay, which are related by the transfer function ip(t), the measured emission-line flux 
can be written as 

/oo 
i>(r)C(t-T)dT (Al) 
-oo 

This equation can be also written using the CCF(t), ACF cor j t (t) and the transfer function 
ip(t) in the form 

/oo 
^(r)ACF cont (t-r)dT (A2) 
-oo 

as has been shown by Penston (1991), Koratkar & Gaskell (1991), and Peterson (1993, 2000). 
Hence, the measured time delay depends on both the ACF cont of the AGN continuum and 
the transfer function ip(t). While the width of the ACF cont at full-width-at-half maximum 
(FWHM) of the continuum variations studied by Sergeev et al. (2002,2011) and Shapovalova 
et al. (2010) amounts to about FWHM~ 1600 days to 2000 days, the FWHM of the contin- 
uum ACF in our 2005 monitoring campaign is FWHM(ACF cont ) ~ 27 days, i.e., about ~ 60 
to ~ 75 times narrower. In the context of investigating the impact of variations on longer 
time scales and the necessity of detrending light curves for those trends to recover the time 
delay, Welsh (1999) mentioned that the points in the ACF and CCF are highly correlated 
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(e.g., Jenkins & Watts 1969) and that those correlations result in spurious large values of 
the CCF, especially for light curves which are characterized by intrinsically broad peaks. We 
suspect that the different time delays are predominantly caused by the significantly different 
widths of the ACF cont (cf., §5.1). 

In addition, using the variations of the broad emission-line flux over times scales of the 
order of the dynamical time will include not only variations which are directly associated with 
continuum variations, i.e., pure reverberation events, but will also include variations caused 
by changes in the distribution and conditions of the line emitting gas. Those variations 
are manifested, for example, in emission-line profile changes which are uncorrelated with 
continuum variability, as shown by Wanders & Peterson (1996). Hence, including emission- 
line profile variability will dilute a reverberation signal, i.e., the direct response of the gas 
on continuum variations. This aspect needs to be addressed when studying the relation of 
continuum and emission-line variability over time scales which are comparable or even longer 
than the dynamical time scale for an AGN. 

Based on the variability of NGC 5548 which has currently the best covered spectroscop- 
ically variability history for the optical continuum and the broad H/3 emission-line flux, it 
has been shown that the measured time delay is strongly correlated with the strength of 
the continuum, i.e., whether the AGN is in a high or low state. Depending on the contin- 
uum luminosity, the delay of the H/3 flux response to continuum variations varies between 
t(H/3) = 4.2 days and t(H/3) = 26.4 days (Bentz et al. 2009b; Peterson et al.2002). The 
continuum strength in the 13 years which were investigated in these studies varied by a 
factor of about ~ 10. The light curves which were studied by Sergeev et al. (2002, 2011) and 
Shapovalova et al. (2010) cover more than 10 years. During this time the strength of the 
continuum of 3C 390.3 varied by a factor of nearly ~ 6. 

Finally, Home et al. (2004) have studied the impact of the duration of a monitoring 
campaign on the measured time delay. They provide a guide line for the optimal length of 
monitoring campaign depending on the brightness of an AGN, to optimize also the use of 
telescope time. They found that a campaign should last at least about 3 times the light 
crossing time of the BLR to recover the velocity resolved transfer function i/j(t : v). In the 
case of cross-correlation functions, shorter campaigns can still yield reliable time delays, in 
particular if the light curve displays features of increasing and decreasing continuum and 
emission-line flux and not only a monotonic increase or decrease of the continuum and 
emission-line flux strength. 
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Figure 1 - A typical g-band image of the field of 3C 390.3, taken with RetroCam on 
September 25, 2005 (ti nt = 120 sec). 3C 390.3 and the comparison stars are marked. In 
addition, the aperture and the location and size of the region which we used to correct for 

the sky background, is shown. 
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Figure 2 - Decomposition of a spectrum of 3C 390.3, observed on November 11, 2005. In 
the top panel the rest-frame spectrum is shown together with the power law continuum fit 
(dashed line), the host galaxy fit (upper thin solid line), the optical Fell emission fit (lower 
thin solid line), and the resulting combined fit (red line). In the bottom panel the residual 
spectrum is displayed which shows the pure emission-line spectrum of 3C 390.3. 
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Figure 3 - The spectrum of 3C 390.3 during a low state (August 30, 1994, 2.2 m at Calar 

Alto Observatory/ Spain). 
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Figure 4 - The H/3 - [O m]AA4959, 5007 emission- line complex of the spectrum shown in 
Fig. 3 (thick line). The scaled Ha profile is shown as a long-dashed line. The residual 
spectrum, which is used to fit the [O ill] emission-line profiles, is shown as the red line. 



-43- 




Figure 5 - The fit of the [O ill] A4959 and [O in] A5007 emission-line profiles employing a 
strong narrow and a weaker, slightly blue-shifted broader Gaussian component for each line 
(thin black lines). The resulting fit is shown as a thick red line which recovers the observed 
spectrum well, as can be seen by the residuum (dashed line). 
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Figure 6 - Light curves of 3C 390.3 for the <^band (in g-band magnitudes), the Fa(5100 A) 
AGN continuum (in units of 10~ 15 ergs -1 cm~ 2 A -1 ) and the variations of the broad 
emission-line fluxes of Ha, H/3, H7, and HellA4686 (displayed in units of 
^Q-i5 er g g -i cm -2^ t]-^ uncertainties of the g-band magnitudes are smaller than the plot 

symbol. 
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Figure 7 - In the top panel the H/3 - [O lll]AA4959, 5007 complex is shown together with 
the scaled broad, double-peaked profile fit of Ha for the H/3 and H7 emission lines (thick 
line). In addition, the profile fits of the narrow emission lines for this wavelength range are 
displayed. In the bottom panel, the residual spectrum is presented after subtracting the 
profile fits of the narrow emission lines and the broad profile fits for H/3 and H7. The 
residual profile at the location of the He II A4686 emission line is reconstructed with a 
scaled double-peaked profile (dotted line) and with a Gaussian profile (thick line). 
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Figure 8 - The Gaussian profile fit for the HellA4686 emission-line profile, under the 
assumption that the blue-wing can be used to represent the red wing (thick solid lines). 
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Figure 9 - The mean and rms spectra of the hydrogen Balmer lines Ha, H/3, and H7. The 
top panel shows the mean spectra (thick lines) and the narrow-line subtracted mean 
spectra (thin lines). In the bottom panel the rms spectrum is presented (solid line) 
together with the location of narrow emission lines (dotted lines). 
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Figure 10 - The mean spectrum of the H/3 4861 emission line is shown (top panel), 
uncorrected and corrected for contamination by He IIA4686 emission. In addition, the mean 
narrow and broad HellA4686 line profile is plotted (thin solid line). In the bottom panel 

the corresponding rms spectra of the H/3A4861 line are shown (thick solid lines). The 
difference between these rms spectra can be associated with the He IIA4686 rms spectrum 
(thin line). The Gaussian profile shaped rms spectrum for HellA4686 is shown for 

comparison (thin solid line). 



-49 - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




= I i i i I i i i I i i i I i i i I L 



4400 4600 4800 5000 5200 

wav elength [A] 

Figure 11 - The H/3 A4861 rms spectrum (thin solid line) and the re-binned, as well as the 
re- binned and scaled Ha A6563 rms spectrum (thick solid line) are shown (top panel). In 
the bottom panel the difference of the H/3 and the re-binned and scaled Ha rms spectra is 

displayed (thick solid line). For comparison, the HeIlA4686 rms spectrum as shown in 
Fig. 10 is shown as thin solid line, as well as the Gaussian profile shaped rms spectrum for 

HellA4686 (dashed line). 
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Figure 12 - The cross-correlation function of the g-h&nd flux variations with the variable 

AGN continuum F c (5100 A) as the driving light curve. The location of the centroid 
(r > 0.8 CCF max ) is shown as dashed line, while the peak of the ICCF is marked by the 

dashed-dotted line. 
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Figure 13 - The cross-correlation functions of the broad emission lines of Ha, H/3, H7, and 
He IIA4686 (thick line) with the g-band variations as the driving light curve. For 
comparison, the thin lines show the cross-correlation functions using the continuum 
variations at A = 5100 A as driving light curve. The locations of the centroid for a 
threshold of 30% 50% and 80% are shown, respectively, for the ICCF based on the <^band 

variation (filled diamonds). 
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Figure 14 - The location of the extraction windows for the blue and red wing, the blue and 
red peak, and the center of the broad line profiles of the Ha (top panel) and H/3 (bottom 
panel) emission lines using the mean spectrum (solid line) and the rms spectrum (dotted 
line, scaled by a factor of 20). The long dashed line indicates the measured location of the 

blue and red peak. 
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Figure 15 - The cross-correlation functions using the AGN continuum Fa (5 100 A) as driver 

light curve (solid line) for the blue wing, blue peak, line center, red peak, and red wing 
(top to bottom) of the broad Ha line profile (on the left) and of the broad H/3 line profile 
(on the right). The ICCFs obtained with the g-band variations as driving continuum are 

shown as dashed line. The location of the centroid, based on a threshold of 80 % of 
lCCF max and the corresponding uncertainty is shown as filled diamonds (F A (5100A)) and 

open diamonds ((/-band), respectively. 
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Figure 16 - Cross-correlation functions for Ha (left panels) and for H/3 (right panels) 
correlating the blue and red wing, as well as the blue and red peak of the line profiles 
directly, with the blue light curve as driver and the red light curve in response (top two 
panels). The following panels show the ICCFs using the variations of the line center as 
driving light curve and the profile wings and peaks in response. 
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Figure 17 - Comparison of the black hole mass estimates that are based on the measured 
time delay for the Balmer lines and the helium He II A4686 line using the second moment 
(June of the line profiles (top panel) and the blue - red peak separation (bottom panel). The 
black hole masses based on the mean line profiles and of the blue - red peak separations 
are shown as filled symbols while the rms spectra results are given by the open symbols. 
The diamonds represent the black hole masses using the AGN continuum F A (5100A) to 
determine the time delay while the boxes display the results using the <^band variations. 
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Figure 18 - Comparison of the normalized light curves of the AGN continuum and of the 
broad H/3 emission line flux of this study (black boxes), of Sergeev et al. (2002,2011, blue 
diamonds), and of Shapovalova et al. (2010, red triangles). The continuum is shown as 
filled symbols and the H/3 emission-line flux is shown as open symbols. For better 
comparison the normalized light curves have been offset. 
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Table 1. Imaging and Spectroscopy Observing Log of 3C 390.3. 



civil date Julian Date - 245 0000 
mm-dd-yy photom. spectros. 
(1) (2) (3) 



09 


16 


05 


3630 


.71895 


3630 


.74012 


09 


17 


05 


3631 


.68545 


3631 


.71968 


09 


18 


05 


3632 


.67934 


3632 


.69208 


09 


22 


05 






3636 


.67738 


09 


24 


05 


3638 


.60241 


3638 


.61469 


09 


25 


05 


3639 


.60314 


3639 


.61829 


09 


30 


05 


3644 


.61169 


3644 


.63453 


10 


02 


05 


3646 


.58313 


3646 


.60055 


10 


10 


05 


3654 


.65719 


3654 


.64235 


10 


11 


05 


3655 


.60982 


3655 


.59675 


10 


12 


05 


3656 


.60524 


3656 


.59131 


10 


17 


05 






3661 


.61921 


10 


19 


05 


3663 


.61128 


3663 


.59896 


10 


20 


05 


3664 


.62103 


3664 


.60890 


10 


21 


05 


3665 


.60079 


3665 


.58876 


10 


23 


05 


3667 


.56875 


3667 


.58457 


11 


01 


05 


3676 


.57946 


3676 


.59541 


11 


03 


05 


3678 


.59495 


3678 


.58313 


11 


11 


05 


3686 


.58735 


3686 


.57751 


11 


12 


05 


3687 


.57943 


3687 


.56689 


11 


13 


05 


3688 


.57791 


3688 


.56608 


11 


21 


05 


3696 


.55431 


3696 


.57008 


11 


22 


05 






3697 


.63839 


11 


24 


05 


3699 


.61123 


3699 


.63440 


11 


25 


05 


3700 


.55593 


3700 


.56501 


11 


27 


05 


3702 


.57947 


3702 


.56630 


12 


06 


05 


3711 


.59450 


3711 


.58072 


12 


08 


05 


3713 


.56451 


3713 


.58304 



a beginning of the night 



-59 - 



Table 2. Mean NLR emission-line rest-frame fluxes, relative to the flux of H/3. 



line flux flux cor 



(1) 




( 


2) 








(3) 






[On]A3727 


1. 


18 


± 


0. 


08 


1. 


41 


± 


0. 


10 


[Nem]A3869 


1. 


20 


± 


0. 


.08 


1. 


41 


± 


0. 


10 


[Nem]A3968 


0. 


.73 


± 





.02 


0. 


.84 


± 


0. 


02 


H7 A4340 


0. 


44 


± 


0. 


02 


0. 


.48 


± 


0. 


02 


[O in]A4363 


0. 


.85 


± 


0. 


02 


0. 


92 


± 


0. 


03 


HcnA4686 


0. 


15 


± 





.01 


0. 


15 


± 


0. 


01 


H/3 A4861 


1. 


00 


± 


0. 


01 


1. 


00 


± 


0. 


02 


[Oni]A4959 


2. 


.48 


± 


0. 


02 


2. 


.44 


± 


0. 


05 


[Om]A5007 


7. 


30 


± 


0. 


.07 


7. 


.12 


± 


0. 


15 


[Fcvn]A5721 


0. 


15 


± 


0. 


.01 


0. 


13 


± 


0. 


01 


He IA5876 


0. 


15 


± 


0. 


01 


0. 


13 


± 


0. 


01 


[Fcvn]A6087 





36 


± 





.02 


0. 


30 


± 


0. 


02 


[Oi]A6300 


0. 


99 


± 





.01 


0. 


.82 


± 





02 


[O i]A6364 


0. 


33 


± 


0. 


01 


0. 


27 


± 


0. 


01 


[Fcx]A6374 


0. 


06 


± 


0. 


.01 


0. 


05 


± 


0. 


01 


[Nn]A6548 


0. 


44 


± 


0. 


.01 





35 


± 


0. 


01 


Ha A6563 


3. 


53 


± 


0. 


07 


2. 


85 


± 


0. 


08 


[Nn]A6583 


1. 


34 


± 


0. 


.04 


1. 


08 


± 


0. 


04 


[Sn]A6716 


0. 


20 


± 





.01 


0. 


16 


± 


0. 


01 


[Sn]A6731 


0. 


22 


± 


0. 


01 


0. 


18 


± 


0. 


01 



integrated H/3-narrow line flux, F(H/3) = 
(26.50 ±0.19) X 10 -15 crgs-^m" 2 

b line ratios are reddening corrected as- 
suming a Balmer decrement of F(Hct) /F(H/3) 
= 2.85. 



Table 3. Observed g-hand magnitudes, measured rest-frame continuum flux density at 
A = 5100 A, and integrated flux of the broad emission lines of Ha, H/3, H7, and HellA4686. 



Julian Date g-band Julian Date F A 5100A a Ha A6563 b H/3A4861 6 H7A434O 6 HeIlA4686' 

-245 0000 [mag] -245 0000 

(1) (2) (3) (4) (5) (6) (7) (8) 



3630 


71895 


14 


543 


± 


001 


3630 


74012 


6 


66 ± 


07 


2123 


4 


± 


33 


7 


617.7 ± 13.0 


114 


2 


± 


17 





HI r* { 1 p" T" 

91.6 ± 15.7 


3631 


68545 


14 


559 


± 


001 


3631 


71968 


6 


66 ± 


07 


2203 


3 


± 


10 


3 


566.6 ± 17.8 


125 


5 


± 


17 


1 


85.0 ± 18.8 


3632 


67934 


14 


566 


± 


001 


3632 


69208 


6 


63 ± 


07 


2457 


7 


± 


35 


3 


p" p" t p" 1 1 r\ 

557.5 ± 12.0 


135 


8 


± 


13 


1 


65.9 ± 13.7 














3636 


67738 


6 


63 ± 


07 


2433 


4 


± 


45 


7 


tr on 1 ac a 

580.2 ± 16.4 


135 





± 


17 


2 


A O A _1_ 1 T *7 

48.4 ± 17.7 


3638 


60241 


14 


589 


± 


001 


3638 


61469 


6 


60 ± 


07 


2379 


1 


± 


43 


4 


t* r' p" 1 1 c\ c\ 

562.5 ± 10.9 


128 


2 


± 


17 


1 


59.1 ± 12.4 


3639 


60314 


14 


590 


± 


001 


3639 


61829 


6 


60 ± 


07 


2320 


6 


± 


29 


5 


P" O A A 1 "1 O "1 

584.4 ± 13.1 


130 


1 


± 


17 


1 


TO A 1 1 a r 

78.4 ± 14.5 


3644 


61169 


14 


598 


± 


001 


3644 


63453 


6 


54 ± 


07 


2368 


3 


± 


29 


3 


586.6 ± 11.4 


127 


8 


± 


14 





80.0 ± 13.3 


3646 


58313 


14 


583 


± 


001 


3646 


60055 


6 


70 ± 


07 


2642 


3 


± 


39 


7 


616.4 ± 11.3 


125 


4 


± 


15 


9 


89.1 ± 14.0 


3654 


65719 


14 


566 


±0 


001 


3654 


64235 


6 


71 ±0 


07 


2545 


6 


± 


38 


1 


583.7 ± 10.2 


127 





± 


17 


2 


77.6 ± 12.0 


3655 


60982 


14 


559 


±0 


001 


3655 


59675 


6 


92 ±0 


07 


2410 


9 


± 


33 


3 


538.2 ±9.8 


133 


5 


± 


15 


3 


67.2 ± 10.5 


3656 


60524 


14 


549 


±0 


001 


3656 


59131 


7 


01 ±0 


07 


2437 


8 


± 


47 


5 


576.2 ± 19.5 


131 


9 


± 


17 


3 


66.3 ±20.4 














3661 


61921 


6 


99 ±0 


07 


2673 


7 


± 


49 


6 


585.2 ±27.9 


150 


2 


± 


35 


8 


51.0 ±28.7 


3663 


61128 


14 


526 


±0 


001 


3663 


59896 


7 


14 ±0 


07 


2627 


9 


± 


47 


4 


572.5 ± 11.3 


159 


4 


± 


12 


8 


67.2 ± 12.1 


3664 


62103 


14 


525 


±0 


002 


3664 


60890 


7 


27 ±0 


07 


2488 


4 


± 


33 


6 


570.2 ± 12.4 


156 


6 


± 


16 


6 


73.9 ± 13.0 


3665 


60079 


14 


522 


±0 


002 


3665 


58876 


7 


26 ±0 


08 


2534 


5 


± 


34 


4 


574.9 ± 15.2 


150 


2 


± 


16 


6 


79.8 ± 15.7 


3667 


56875 


14 


495 


±0 


004 


3667 


58457 


7 


31 ±0 


08 


2611 


2 


± 


52 


2 


598.1 ±27.4 


151 


4 


± 


15 


5 


84.0 ±28.0 


3676 


57946 


14 


520 


±0 


002 


3676 


59541 


7 


20 ±0 


08 


2460 


7 


± 


46 


3 


605.6 ±9.2 


160 


6 


± 


16 


9 


99.0 ± 10.7 


3678 


59495 


14 


529 


±0 


002 


3678 


58313 


6 


94 ±0 


07 


2488 


4 


± 


37 


3 


623.8 ± 14.0 


152 


5 


± 


17 


5 


116.7 ± 15.9 


3686 


58735 


14 


574 


±0 


001 


3686 


57751 


6 


78 ±0 


07 


2334 


8 


± 


40 


1 


610.1 ±6.2 


169 


2 


± 


16 


8 


109.6 ± 7.9 


3687 


57943 


14 


591 


±0 


002 


3687 


56689 


6 


74 ±0 


07 


2340 


7 


± 


36 


4 


594.4 ± 10.0 


170 


6 


± 


16 


7 


107.7 ± 11.0 


3688 


57791 


14 


592 


±0 


002 


3688 


56608 


6 


70 ±0 


07 


2523 


8 


± 


34 


1 


609.5 ±22.9 


170 


2 


± 


17 





108.3 ±23.7 


3696 


55431 


14 


633 


±0 


001 


3696 


57008 


6 


24 ±0 


07 


2291 


8 


± 


30 





619.6 ± 13.3 


169 


8 


± 


16 


9 


110.2 ± 14.3 














3697 


63839 


6 


21 ±0 


07 


2162 


3 


± 


40 


4 


632.8 ± 15.5 


145 


3 


± 


20 


6 


110.2 ± 15.6 


3699 


61123 


14 


651 


±0 


002 


3699 


63440 


6 


09 ±0 


06 


2427 


7 


± 


33 


6 


624.6 ±8.6 


170 


5 


± 


17 


1 


123.3 ± 10.6 


3700 


55593 


14 


647 


±0 


002 


3700 


56501 


5 


89 ±0 


06 


2399 


6 


± 


40 


6 


634.3 ±21.1 


165 


6 


± 


17 


4 


113.0 ±22.3 


3702 


57947 


14 


659 


±0 


002 


3702 


56630 


5 


91 ±0 


06 


2372 


6 


± 


36 


2 


621.5 ± 12.2 


162 


2 


± 


17 


2 


111.2 ± 13.8 


3711 


59450 


14 


708 


±0 


002 


3711 


58072 


5 


65 ±0 


06 


2441 


6 


± 


47 


4 


676.9 ± 15.0 


171 





± 


12 


8 


104.8 ± 16.2 


3713 


56451 


14 


694 


±0 


002 


3713 


58304 


5 


70 ±0 


06 


2519 


7 


± 


31 


8 


663.0 ± 12.0 


170 


7 


± 


17 


1 


103.9 ± 13.5 



a continuum flux in [10 15 ergs 1 cm 2 A x ] 
b integrated emission-line flux in [10 -15 ergs _1 cm~ 2 ] 
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Table 4. Measurements of the mean and rms line profile parameters of the Balmer lines 

Ha, H/3, H7, and the helium line He IIA4686. 







Mean Spectra 




feature 


FWHM 




FWHM 






[A] 


[A] 


[kms- 1 ] 


[kms" 1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


Ha A 6563 


275 ± 1 


101 ± 1 


12563 ± 31 


4607 ± 29 


H/3A4861 


214 ±1 


87 ±1 


13211 ±28 


5377 ± 37 


H7 A 4340 


172 ±2 


58 ±1 


11875 ± 103 


3991 ± 29 


He 11 A4686 


242 ±8 


103 ±3 


15480 ± 512 


6590 ± 192 






rms Spect 


ra 






FWHM 




FWHM 






[A] 


[A] 


[kms- 1 ] 


[km S -1 ] 


Ha A 6563 


278 ± 33 


106 ±6 


12679 ± 1312 


4839 ± 215 


H/3A4861 


176 ± 26 


88 ±5 


10872 ± 1670 


5455 ± 278 


H 7 A4340 


165 ± 19 


75 ±2 


11422 ± 1458 


5191 ±82 


He 11 A4686 


207 ± 15 


80 ±11 


13244 ± 960 


5125 ± 704 



Table 5. Statistical properties of the continuum and broad emission-line variations 



feature mean flux c 
(1) (2) 


<4 
(3) 


Fuar 

(4) 


min 

(5) 


max 
(6) 


g-band a 


14.58 


0.28 


0.050 


14.71 


14.50 


V-band b 


7.34 


0.37 


0.051 


6.53 


7.94 


F A (5100A)<> 


6.63 


0.47 


0.070 


5.65 


7.31 


Ha A 6563 


2439 


128 


0.041 


2123 


2674 


H/3 A 4861 


598 


32 


0.047 


538 


677 


H7 A 4340 


149 


18 


0.031 


114 


171 


He 11 A4686 


88 


21 


0.155 


48 


123 


a continuum flux 


in [mag] 










b continuum flux 


in [10- 15 


ergs -1 


cm -2 A 






c integrated line 


flux in [10 


-is cr g 


s _1 cm - 


2 ] 
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Table 6. Results of the cross-correlation analysis of the broad line emission of the Balmer 

lines Ha, H/3, H7, the helium line HellA4686, using AGN continuum variations at 
F A (5100A) and the <^band variations as driving continuum. The centroid of the ICCF has 
been computed for 80 % of the corresponding maximum of the ICCF. The delays calculated 

with SPEAR are given in column (8). 



feature 




F A (5100A) 






g- band 




SPEAR 




Tcent 


Tpeak 


CCF max 


7~cent 


T peak 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



g-band 


o-2±i:i 




0.97 ±0.01 








55+ 1 - 56 

U.OO_ 45 


Ha A 6563 


56.3±*;g 




0.83 ± 0.39 


56.3l?:| 


56.5l 2 6 :° 


0.87 ±0.30 


JZ - J -0.6 


H/3A4861 


46.4j£° 




0.91 ± 0.07 


44.31*;° 


45.0^-0 


0.87 ±0.10 




H7A4340 


01 7 +11.8 


32-0±^J 


0.75 ±0.11 


58.ll|:? 


58.0l« 


0.86 ± 0.34 


32.1 ± 17.3 


HeiiA4686 


24.5l^ 7 


24 5+ 6 ' 


0.75 ± 0.11 


22.3tH 


23-011° 


0.67 ±0.13 


36.0 ± 5.2 



Table 7. Location of the blue and red peak, as well as of a red feature in the broad line 
profiles of the Ha, H/3, and H7 emission, if detectable. 









Mean Spectra 








feature 


blue peak 


red peak 


red feature 




[A] 


[kms- 1 ] 


[A] 


[kms- 1 ] 


[A] 


[kms- 1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Ha A 6563 


6488.3 ± 1.0 


-3413 ± 41 


6662.1 ±0.0 


±4538 ± 54 


6613.0 ±1.0 


+2292 ± 38 


H/3A4861 


4807.4 ±1.0 


-3324 ± 49 


4935.3 ± 0.0 


±4563 ± 68 


4901.6 ± 1.0 


+2485 ± 42 


H7A4340 


4287.6 ± 1.0 


-3650 ± 69 


4409.7 ± 1.5 


±4783 ± 104 












rms Spectra 








feature 


blue peak 


red peak 


red feature 




[A] 


[kms- 1 ] 


[A] 


[kms- 1 ] 


[A] 


[kms- 1 ] 


Ha A 6563 


6483.5 ± 2.0 


-3622 ± 91 


6665.9 ± 2.0 


+4711 ± 228 


6619.3 ±4.0 


+2582 ± 183 


H/3A4861 


4801.5 ±3.0 


-3689 ± 185 


4945.8 ± 4.0 


+5210 ± 247 


4897.1 ±6.0 


±2207 ± 370 


H 7 A4340 


4293.0 ± 10.0 


-3277 ±690 


4402.0 ± 6.0 


+4251 ± 415 
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Table 8. Range of the line profile regions to extract light curves for the profile wings, the 
line center, and the blue and read peak in the broad emission- line profiles of the Ha and 
H/3 lines. In brackets the velocity-ranges are given which are used by Gezari et al. (2007). 



feature 
(1) 


Ha A 6563 
(2) 




H/3A4861 
(3) 




blue wing 


6345.5 A— 


6455.0 A 




4702.0 A- 4783.0 A 






-9916 kms" 1 - 


-4916 km s 


-l 


-9819 km s" 1 4819 km s 


-1 




(-12300 km s- 1 - 


- -4500 km £ 








blue peak 


6455.0 A— 


6521.0 A 




4783.0 A- 4832.0 A 






-4916 km s" 1 - 


-1916 kins 


-l 


-4819 km s" 1 1819 km s 


-1 




(-4500 km s" 1 - 


-1700 km s 


- 1 ) 






center 


6521.0 A— 


6629.0 A 




4832.0 A- 4911.0 A 






-mekms- 1 - 


+3032 km s 


-l 


-lSigkms- 1 f-3063kms 


-1 




(-1700 km s" 1 - 


+3300 km s 


- 1 ) 






red peak 


6629.0 A— 


6695.0 A 




4911. OA- 4959.5A 






+3032 km s" 1 - 


+6032 km s 


-l 


+3063 kins" 1 h6063kms 


-1 




(+3300 km s- 1 - 


+5900 km s 


- 1 ) 






red wing 


6695.0 A— 


6804.5 A 




4959.5 A- 5040.5 A 






+6032 km s" 1 - 


+ 11032 km e 


-i 


+6063 km s" 1 hi 1063 km 


s- 1 




(+5900 km s- 1 - 


+15000 km s 


i- 1 ) 
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Table 9. Results of the cross-correlation analysis of the wings, the center, and the blue 
and red peak of the broad Balmer lines Ha and H/3 using AGN continuum variations at 
Fa(5100A) and the g-band variations as driving continuum. The centroid of the ICCF has 
been computed for 80 % of the corresponding maximum of the ICCF. 



feature 




F A (5100A) 






#-band 






Tcent 


Tpeak 


CCF max 


Tcent 


T peak 


CCFmax 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



Ha 


blue wing 




58.0l^ 





88 ± 


.47 


55.91*;° 


57.0l 3 5 ;° 


0. 


92 


± 





59 


Ha 


blue peak 


54.3±2"2 


56.5l« 





84 ± 


.49 


53.2l«;J 


55.5 to 





85 


± 





56 


Ha 


center 


54.5t 6 4 « 


56.5l« 





82 ± 0. 


.41 


53.5llg 


56. Ol" 





85 


± 





73 


Ha 


red peak 


54.lt 6 ,l 


56.5l« 


0. 


.85 ± 0. 


26 


52.7tU 


55.01*'° 





86 


± 


0. 


20 


Ha 


red wing 




57.511:° 





89 ± 0. 


.40 


49.91JV 


51-0li° 5 


0. 


.95 


± 


1 


.28 


Ha 


blue vs. red wing 


-3-415:1 


-2 5+ 7 - 8 


0. 


.64 ± 0. 


10 
















Ha 


blue vs. red peak 




0+1.0 
u - u -0.5 


0. 


.67 ±0. 


10 

















H/3 


blue wing 


49.611?* 


52 n+ 100 





.84 ± 


.34 


48.4l^ 8 3 


49.5lJJ;° 





.84 


±0 


16 


H/3 


blue peak 


48.0+^ 3 7 3 


oi - u -14.5 





.86 ± 


71 


47.2+J^ 1 


49.olii;° 





85 


±0 


.27 


H/3 


center 


49.6l 1 1 2 3 : 4 r 


53 n+ 90 


0. 


.85 ± 0. 


56 


49-OlJol 


50 5+ 10 5 
ou.o_ 12 


0. 


84 


±0 


19 


H/3 


red peak 


5i.6i? : g 


54.5l« 


0. 


.80 ± 0. 


19 


53.11^ 


■54 5+6-0 


0. 


86 


±0 


74 


H/3 


red wing 


54.4l? °3 


56.01" 


0. 


.81 ±0. 


35 


56.3l|l 


57.0l 5 8 ;5 


0. 


91 


±0 


71 


H/3 


blue vs. red wing 


2.9±ll:l 


o n +i8.o 

• 3 - u -22.0 





47 ± 0. 


10 














H/3 


blue vs. red peak 


i-ol|:I 


0+ 100 
u - u -io.o 





.50 ±0 


11 
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Table 10. Estimates of the black hole virial product M^ r , based on the Balmer and 
helium line profile properties [oime and separation of the blue and red peak in the 
double-peaked profiles) and the time delays r cent (80 % threshold) of the broad 

emission-line flux variations. 



Virial Black Hole Product M 



bh 



in 10 s Mr 



feature 
(1) 



using au ne of the line profiles 

F A (5100A) 
M" r (mean) 



(2) 



M" r (rms) 



(3) 



M™ r (mean) 



(j-band 



('I) 



M" r (rms) 



(5) 



Ha A 6563 


9 OO+0.17 
z -°°-0.16 


2.57" 


H/3A4861 


9 fi9 +0.20 
-0.18 


2.69; 


H7 A 4340 




i.67; 


HcnA4686 


2.08l°;« 


1.26 


average 


2-oitS:S 


2.05; 



+0.25 
-0.24 
+0.29 
-0.27 
+0.62 
-0.55 
+0.36 
-0.35 

+0.38 
-0.36 



2.33 
2.50 
1.81 
1.89 



0.10 
0.27 
+0.17 
0.19 
+0.14 
0.19 
0.56 
-0.33 



2.57 

2.57 
3.05 
1.14 



0.20 
0.34 
+0.25 
0.27 
+0.23 
0.33 
+0.40 
-0.30 



o in+0.24 
Z - ±J -0.21 



z.oo_ 37 



using the separation of the blue and red peak 
feature F A (5100A) g-band 



(1) 



(2) 



(3) 



(4) 



(5) 



Ha A 6563 
H/3A4861 
H7 A 4340 
average 



1.74 
1.41 
1.10 
1.42 



+0.14 
0.13 
+0.12 
0.10 
+0.41 
-0.37 

+0.35 
0.32 



1 qi+0-21 
i - yl -0.21 
+0.23 
0.22 
+0.42 
-0.39 



1.87 
0.88 



1.55 



+0.48 



1.74 
1.34 
2.02 

1.70 



+0.09 
0.21 
+0.10 
0.11 
+0.19 
-0.23 

-0.26 
-0.33 



1 qi+0.18 

1 7 o+0.21 
1 -'°-0.22 



1.61 



+0.50 
0.51 



1.77 



+0.29 
-0.31 



